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NOMENCLATURE 


Symbol  Description  Units 


A 

cross  sectional  area  of  aerosol  conduit 

ft^ 

A 

c 

total  cross  sectional  area  of  particles 
representing  a  size  distribution 

ft^ 

(A  )  ^ 

^  c  nd 

cross  sectional  area  of  n  particles  with 
diameter  d 

ft2 

A 

s 

total  surface  area  of  particles  representing 
a  size  distribution 

ft^ 

A 

w 

wall  area  of  aerosol  conduit 

ft^ 

a 

subscript  meaning  air 

- 

a 

particle  radius 

a 

constant  for  X 

- 

atm 

atmosphere 

- 

b 

constant  for  X 

- 

c 

constant  for  X 

- 

C 

particle  concentration  by  volume,  (ft  solids/ 

ft^  of  aerosol 

dimensionless 

C 

n 

particle  number  concentration 

C 

average  air  heat  capacity 

Btu/lb 

°F 

CFM 

cubic  feet  per  minute 

- 

C 

P 

particle  heat  capacity 

Btu/lb 

°F 

c 

p 

average  particle  heat  capacity 

Btu/lb 

°F 

dA 

differential  area  of  conduit  wall 

positioned  with  respect  to  dV 

ft 

V 


d 

P 

particle  diameter 

ft 

dV 

differential  volume  of  aerosol  arbitrarily 
positioned  and  receiving  radiation  from  dA^ 

ft^ 

ordinate  for  figures  32  through  37;  heat 
absorbed  per  unit  wall  radiation  intensity 
per  unit  cross  sectional  area  of  enclosure 

dimensionless 

H 

enthalpy  per  unit  time 

Btu/min 

h 

a 

convection  heat  transfer  coefficient  for  air, 
and  from  a  particle  to  its  surroundings 

Btu/hr  ft^  °F 

h 

P 

convection  heat  transfer  coefficient  for  a 
particle  suspension 

Btu/hr  ft^  °F 

I 

denotes  integral  of  equation  (2.19) 

dimensionless 

I 

O' 

alpha  series  of  I 

dimensionless 

^al 

first  term  of  alpha  series 

dimensionless 

second  term  of  alpha  series 

dimensionless 

I  o 
<x3 

all  odd  terms  of  alpha  series  for  n  >  1 

dimensionless 

^ah 

all  even  terms  of  alpha  series  for  n  >  2 

dimensionless 

h 

beta  series  of  I 

dimensionless 

hi 

first  term  of  beta  series 

dimensionless 

h2 

second  term  of  beta  series 

dimensionless 

h,n  >  2 

all  terms  of  beta  series  for  n  >  2 

dimensionless 

i 

summation  index 

- 

black  body  radiation  intensity  normal  to 
an  emitting  surface 

Btu/hr  ft^  (ju 

^bnl 

black  body  radiation  intensity  normal  to 
surface  1 

Btu/hr  ft^  (JO 

\n2 

black  body  radiation  intensity  normal  to 
surface  2 

Btu/hr  ft^  (JO 

^bnw 

black  body  radiation  intensity  emitted 
radially  inward  from  a  cylindrical  wall 

Btu/hr  ft^  (JO 

vi 

i 

o 


iDlack  body  radiation  emitted  from  a  cylinder 

wall  at  an  angle  with  the  radial  direction  at  2 

the  point  of  emission  Btu/hr  ft  co 

directional  radiation  intensity  at  reference  2 

point  Btu/hr  ft  u) 


i 

X 


radiation  intensity  at  point  x  in  the  x 
direction 


Btu/hr  ft 


2 


UD 


J 


integral  term  defined  by  equation  (2.43)  dimensionless 


J 


1 


J-terms  not  having  a  k  coefficient; 
equation  (2.43). 


J 


E 


sum  of  J-terms  having  a  k  coefficient; 
equation  (2.43) 


dimensionless 


dimensionless 


J 

k 

k 


k 

P 


L 


m 


n 


P 

q 

P 


suimnation  index 

1*39 

constant  for  the  relation  k/1 
attenuation  coefficient 


attenuation  coefficient  for  gases ^  negative 
number 

thermal  conductivity  of  particle 
thermal  conductivity  of  air 


Btu/hr  ft 
Btu/hr  ft  °F 


heated  length  of  radiating  aerosol 
conduit,  and  mean  beam  length  in 
equation  (l.3) 

geometric  variable  defined  by  equation  (2.10) 
and  described  in  Figure  6. 

general  term  number  of  series  expansion,  number 
of  particles  having  same  diameter,  and  number 
of  moles  of  air  in  aerosol  conduit 


exponent 

isotropic  radiation  flux  Btu/hr  ft 

subscript  indicating  particle 

black  body  radiant  heat  transfer  between  two 

surfaces  without  an  intervening  absorbing 

medium  Btu/min 


vii 


%12 

%21 

% 

%1 


black  body  radiation  emitted  from  surface  1 
and  absorbed  by  surface  2  with  no  intervening 
absorbing  medium 

Btu/min 

black  body  radiation  emitted  from  surface  2 
and  absorbed  by  surface  1  with  no  intervening 
absorbing  medium 

Btu/min 

black  body  energy  absorbed  per  unit  time 
accounting  for  attenuation 

Btu/min 

black  body  energy  absorbed  per  unit  time 
with  attenuation  not  considered 

Btu/min 

black  body  radiation  attenuated  from 

Btu/min 

heat  transferred  by  convection  to  particle 
free  air 

Btu/min 

heat  transferred  by  convection  to  particle 
suspension 

Btu/min 

radius  of  cylinder  and  sphere  and  l/2  distance 
between  parallel  plates 

ft 

variable  radius 

ft 

S  distance  between  an  elemental  area  emitting 

to  a  differential  TOlume  of  aerosol 


ft 


S 

V 


T 

t 

a 

t 

P 

^0 

t 

o 

U 

u 

V 


ratio  of  total  cross  sectional  area  of  particles 
to  the  total  volume  of  particles  composing  a 
complete  distribution 

temperature 

temperature  of  air 

temperature  of  particle  suspension 

temperature  at  time  6 

temperature  at  time  0=0 

overall  heat  transfer  coefficient 

substitution  variable  of  integration 

total  volume  of  particles  composing  a  distribu¬ 
tion,  and  volume  of  heated  aerosol  conduit 


ft 


-1 


Btu/hr  ft^  °F 


ft 


3 


viii 


volume  of  n  particles  of  diameter  d 
substitution  variable  of  integration 


mass  flow  rate 


mass  flow  rate  of  air 


Ibs/min 

Ibs/min 


mass  flow  rate  of  particles 

2 

General  quadratic  function:  X=cx  +  bx  +  a 

length  of  radiating  enclosure  or  length  equal 
to  it;  also  length-to-radius  ratio  for  these 
to  quantities 

geometric  mean  diameter 

axial  distance  separating  differential  aerosol 
volume,  dV,  receiving  radiation  from  emitting 
wall  element,  dAw.  See  Figure  6 

functional  notation  in  Newton-Cotes  formulae 

ordinates  employed  in  Newton-Cotes  formulae 

dummy  variable  of  integration,  equation  (l.4) 

axial  distance  between  differential  aerosol 
volume  and  beginning  of  radiating  aerosol 
conduit 


dimensionless 


GREEK  LETTERS 


particle  absorptivity 

angle  between  direction  of  a  radiant  beam  and 
the  normal  direction  to  the  surface  s 


radians 


angle  between  direction  of  a  radiant  beam 
and  the  normal  direction  to  dAw 


radians 


adsorption  efficiency,  heat  absorbed  divided 
by  total  heat  emitted  i 

difference  operator 

change  in  when  a  length  nx  emits  to  an 
adjacent  length  x.  is  the  converse  for  the 

emitting  and  receiving  length.  ( 


dimensionless 


dimensionless 


ix 


e 

0 

X 

P 

pp 

a 

g 

$ 

Y 


Y 

m 

Y' 

m 


particle  emissivity 
spectral  particle  emissivity 
time 

radiation  wave  length 

unit  of  length 

density  of  air 

particle  density 

standard  geometric  deviation 

p  C  /p 

plane  angle  between  plane  containing  dV  and 
conduit  axis  and  the  plane  containing  dAw 
and  conduit  axis 

value  of  Y  giving  minimum  attenuation 

value  of  Y  as  determined  only  by  Lambert's 
cosine  law 


0)  solid  angle 

Special  notations  for  energy  balance  calculations; 


Superscripts 


description 


0 


no  particle  flow 


1 


particle  flow 


dimensionless 

dimensionless 

sec 

micron 

micron 

Ib/ft^ 

Ib/ft^ 

dimensionless 

dimensionless 

radians 

radians 

radians 

steradian 


Subscripts 

0 

1 

3 

A 

c 

P 


property  evaluated  at  furnace  inlet 
property  evaluated  at  heat  exchanger  inlet 
property  evaluated  at  heat  exchanger  exit 
air 

heat  exchanger  coolant 
particles 
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SUMMARY 


The  principles  describing  the  transfer  of  heat  by  radiation  were 
formulated  in  the  latter  part  of  the  previous  century  and  their  applica¬ 
tion  to  practical  problems  has  been  the  subject  of  much  study  since  that 
time.  The  resulting  analyses  are  mostly  of  limited  applicability,  how¬ 
ever,  because  of  the  lengthy  calculations  involved  and  the  numerous 
empirical  assumptions  required  to  obtain  them.  Wow  space  exploration, 
nuclear  technology,  and  plasma  research,  to  name  a  few,  require  more 
nearly  exact  solutions^  since  the  temperatures  encountered  in  these 
fields  may  be  much  above  those  in  conventional  equipment  and  radiation 
is  likely  the  dominant  mode  of  heat  transfer.  Fortunately,  solutions 
are  attainable  now  that  high-speed  electronic  computers  can  be  employed 
to  solve  problems  that  before  were  hopelessly  involved. 

A  specific  problem  of  timely  concern  is  the  absorption  of  radiant 
energy  by  particles  suspended  in  weakly  absorbing  or  transparent  gases. 
This  investigation  undertook  to  examine  both  theoretically  and  experi¬ 
mentally  the  transference  of  energy  by  radiation  from  a  cylindrical 
enclosure  to  a  particle  cloud  contained  within.  Specifically,  the 
study  employed  a  non-ref lective  wall  radiating  to  essentially  non-re- 
flective  particles.  The  result  was  not  intended  to  be  inclusive;  rather 
it  was  to  be  a  new  approach  using  modern  computational  methods,  which 
could  be  modified  to  describe  successively  more  complex  processes  and 
still  be  reliably  related  to  basic  radiation  principles. 
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The  net  black  body  radiant  interchange  between  two  surfaces  may 


be  described  by 


^^bn2“^bn]9i^  ( 


COS'  3^  cos  p^N 
~2  / 


dA^dA2 


where  i  and  i,  ^  are  the  black  body  radiation  intensities  normal  to 
bnl  bn2 

the  differential  area  dA^  and  dA^,  and  3^  are  the  angles  between 

the  normal  direction  to  the  surfaces  and  the  line  segment  S  connecting 
them.  For  the  particular  problem  investigated  here,  this  expression 
becomes 

2n  2tt  R  L  L-z  ^ 

-J  J  /  J  J 

o  o  o  o  -z 


.2  ^2  2  2 


)  ^^2Rs' 


where  the  radiation  intensity  is  attenuated  exponentially  with  distance 

(e~^)and  the  particles  were  sufficiently  low  in  temperature  not  to 

back  radiate  significantly  to  the  wall,  i.e.,  i^^^  =  0.  The  other  terms 

are:  1,  ,  the  normal  wall  intensity;  L,  the  cylinder  length;  R,  its 

bnw’ 

radius;  y  and  z  (distance);  and  Y  and  0  (angles)  describing  the  geometric 

location  of  a  differential  wall  element  dA  that  radiates  to  a  differen- 

w 

tial  amount  of  aerosol  surface  dA  .  The  integral  heat  transfer  equation 

P 

was  rewritten  as 


R  2n 

%  =  J  J 

o  o 


where 


L  L-z 


dydz 


X  vi 


and  the  0  integration  was  performed  to  give  the  coefficient  2tt.  The 

integral  I  was  evaluated  analytically  hy  means  of  a  series  expansion  of 
Ic  s 

e“  .  The  subsequent  integrations  for  t  and  r  were  performed  numeri¬ 
cally  by  use  of  closed  form  Newton-Cotes  quadrature  formulae.  The 
integral  equation  was  solved  first  for  a  particular  set  of  conditions 
applying  to  a  concurrent  experimental  investigation;  calculations  were 
made  later  to  apply  to  other  systems. 

A  generalized  correlation  was  produced  for  radiation  systems  of 
specific  interest  to  this  investigation  that  related  the  rate  of  heat 
absorption,  radiation  intensity,  cylinder  length  and  radius,  and  parti¬ 
cle  cloud  attenuation  constants.  The  results  are  presented  graphically, 
and  an  interpolation  procedure  is  described  to  obtain  very  nearly  exact 
data  points  not  specifically  represented.  Finally,  equations  were 
derived  for  evaluating  the  radiant  heating  a  particle  cloud  experiences 
while  in  unheated  entrance  and  exhaust  sectiok'^  adjacent  to  a  heated 
length. 

While  the  theoretical  analysis  for  thermal  radiation  processes 
was  treated,  the  results  apply  equally  as  well  to  other  types  of  radia¬ 
tion  such  as  X-rays.  For  example,  the  dose  rate  for  aerosols  in  a 
cylindrical  enclosure  might  be  calculated. 

The  experimental  investigation  sought  to  measure  the  absorption 
of  radiation  by  low  temperature  aerosol  systems,  absorption  data  for 
aerosols  under  these  conditions  not  being  available.  An  experimental 
method  was  developed  so  that  a  particle  cloud  could  be  exposed  to  a 
radiant  field  in  which  conductive  and  convective  heat  transfer  were 
reduced  to  a  minimum  and  the  energy  could  be  absorbed  by  the  cloud  in 
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such  a  way  that  its  magnitude  might  he  determined  by  enthalpy  balances 
before  and  after  exposure.  The  system  exposed  a  particle  cloud  to  a 
radiant  heat  flux  by  passing  it  through  a  quartz  tube  which  transmitted 
to  its  interior  a  large  portion  of  the  radiation  falling  on  it  exter¬ 
nally.  The  quartz  was  protected  from  extraneous  heat  effects  other 
than  radiation  by  enclosing  it  in  an  evacuated  air  space  the  outer 
surface  of  which  was  air  cooled.  The  surrounding  protective  tube  was 
also  quartz.  The  particle  cloud  was  generated  at  a  constant  rate  and 
preconditioned  to  a  constant,  fixed  temperature.  Then  after  passing 
through  the  radiation  field,  the  cloud  was  cooled  calorimetrically  to 
determine  the  heat  absorbed.  The  apparatus  was  designed  for  continuous 
operation  and  a  filter  mechanism  was  developed  to  capture  the  aerosol 
particles  and  return  them  to  the  generator. 

Experimental  data  were  collected  for  aerosols  of  ferrous  sulfide, 
cupric  oxide, and  zinc.  Particle  sizes  were  chiefly  between  20  to  30, 

30  to  44,  44  to  53,  53  to  88  microns  for  the  ferrous  sulfide  and  cupric 
oxide  powders  and  1  to  10  microns  for  the  zinc.  Particle  clouds  of 
these  materials,  when  exposed  to  a  radiant  heat  flux  of  17,600  Btu/hr 
ft  in  a  0,0526  foot  diameter  quartz  conduit  1.0l4  feet  long,  absorbed 
up  to  fifteen  per  cent  of  the  radiation  transmitted  through  the  enclo¬ 
sure  wall. 

The  experimental  data  were  employed  in  testing  the  theoretical 
calculations.  Calculations  of  the  type  employed  had  been  proposed  in 
previous  investigations  and  had  been  applied  in  certain  experimental 
studies  requiring  their  use.  A  direct  test  of  their  validity,  however, 
had  never  before  really  been  performed.  The  experimental  results 
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confirmed  the  theoretical  calculations  and  showed  that  accurate  evalua¬ 
tion  of  the  absorption  of  radiant  energy  by  particles  suspended  within 
a  cylindrical  enclosure  might  be  assessed  if  precise  measurements  of 
the  dispersion  concentration  and  particle  size  and  size  distribution 
data  are  available. 

It  is  recommended  that  subsequent  studies  be  conducted  to  extend 
the  theoretical  development  to  include  back-radiation  and  radiation 
scattering  from  the  particles,  to  investigate  the  effects  of  non-black 
enclosure  surfaces,  and  to  consider  aerosol  temperature  profiles. 
Further  experimental  work  is  recommended  with  systems  having  other 
length-to-radius  ratios  to  test  more  completely  the  theoretical 
development . 
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CHAPTER  I 


INTRODUCTION 

Definition  of  the  Problem 

The  principles  describing  the  transfer  of  heat  hy  radiation  were 
formulated  in  the  latter  part  of  the  previous  century  and  their  applica¬ 
tion  to  practical  problems  has  been  the  subject  of  many  investigations 
since  that  time.  The  majority  of  these,  however,  have  related  to  furnace 
technology  because  here  sufficiently  high  temperatures  were  first  attained 
for  thermal  radiation  processes  to  have  commercial  significance  and  to 
merit  the  scrutiny  of  engineers.  The  problems  were  primarily  associated 
with  the  transfer  of  radiation  between  solid  surfaces  or  between  solid 
surfaces  and  absorbing  gases  such  as  carbon  dioxide  and  water  vapor. 

The  resulting  analyses  were  of  limited  applicability  because  of  the  nu¬ 
merous  empirical  assumptions  about  combustion  rates,  furnace  shapes, 
fuel  characteristics,  stoichiometric  relationships,  etc.,  that  were 
required.  They  were,  nevertheless,  valuable  in  the  design  of  furnaces 
and  aided  immensely  in  evaluating  furnace  performance. 

Radiation  heat  transfer  is  currently  undergoing  a  renewal  of 
interest,  the  advent  of  plasma  generators,  re-entry  missiles,  and  novel 
concepts  in  nuclear  technology  having  brought  to  the  fore  new  and  more 
challenging  radiation  problems.  Solutions  are  needed  to  basic  as  well 
as  practical  problems  in  these  new  areas  where  temperatures  are  much 
above  those  in  conventional  equipment  and  where  radiation  is  the 
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dominant  mode  of  heat  transfer. 

A  specific  problem  of  timely  concern  is  the  absorption  of  radiant 
energy  by  particles  suspended  in  weakly  absorbing  or  transparent  gases. 

In  particular,  one  version  of  a  nuclear  rocket,  as  presently  envisioned, 
requires  the  transfer  of  heat  at  an  extremely  high  rate  to  the  propel¬ 
lant,  which  for  best  performance  should  be  hydrogen  gas.  Since  the  tem¬ 
perature  of  the  source  must  be  exceedingly  high  for  reasons  of  efficiency, 
the  transfer  will  be  primarily  by  radiation.  Many  gases,  including 
hydrogen  below  about  10,000°  F,are  essentially  transparent  however.  It 
has  been  proposed,  therefore,  that  the  gas  be  made  absorptive  below  this 
temperature  by  the  addition  of  small  particles.  This  is  the  immediate 
problem  that  suggested  this  investigation. 

Detailed  information  on  the  absorption  of  radiation  by  aerosols  -- 
the  technical  name  applied  to  a  system  consisting  of  a  gas  with  particles 
dispersed  therein  --  has  applications  in  other  areas  also.  X-ray  and 
nuclear  radiations,  for  instance,  are  known  to  produce  chemical  reactions 
in  dispersed  systems;  just  how  is  of  much  interest  to  the  process  indus¬ 
tries,  for  example.  Particles  can  be  subjected  to  thermal  energy  in 
plasma  generating  devices  to  produce  spherical  particles,  crystal  trans¬ 
formations,  and  size  reductions.  For  protection  in  a  time  of  national 
emergency,  aerosols  might  be  produced  to  absorb  radiation,  both  nuclear 
and  thermal,  that  could  otherwise  prove  lethal. 

Since  high-speed  electronic  computers  can  be  employed  to  solve 
problems  that  before  were  hopelessly  without  solution,  this  investigation 
undertook  to  examine  both  theoretically  and  experimentally  the  specific 
radiation  process  in  which  energy  is  transferred  from  a  cylindrical 


enclosure  to  a  particle  cloud  within.  The  theoretical  portion  sought 
complementary  analytical  solutions  which  were  sufficiently  general  to 
be  employed  in  predicting  the  behavior  of  systems  considerably  differ¬ 
ent  from  the  experimental  one.  The  experimental  study  was  made  to 
produce  direct  basic, data  about  the  process,  unaffected  as  nearly  as 
possible,  by  other  heat  transfer  phenomena.  Specifically,  this  study 
employed  a  non-reflective  wall  radiating  to  essentially  non-reflective 
particles.  The  result  was  not  expected  to  be  inclusive;  rather,  it  was 
to  make  use  of  a  new  approach,  using  modern  computational  methods, 
which  could  be  modified  to  describe  successively  more  complex  processes 
and  still  be  reliably  related  to  basic  radiation  principles. 


Background 

Radiant  heat  transfer  from  a  solid  surface  to  a  cloud  of  parti¬ 
cles  is  similar  to  the  problem  of  radiation  to  an  absorbing  gas.  Since 
the  amount  of  energy  absorbed  by  a  gaseous  medium  is  a  function  of  the 
distance  traversed  by  a  radiant  beam,  the  problems  here  are  the  des¬ 
cription  of  attenuation  processes  in  such  media  and  the  derivation  of 
equations  from  which  to  calculate  the  radiation  properties  of  an 
absorbing  or  emitting  gas  as  a  function  of  its  temperature,  pressure, 
composition,  and  shape.  The  relationship 


(1.1) 


describes  the  process  mathematically  in  most  cases,  or  at  least  can  be 
modified  to  apply  to  them,  i  and  i  are  the  intensities  at  the  surface 
and  at  a  distance  x  from  the  surface.  The  exponential  factor,  k,  is 
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the  absorption  coefficient  which  for  gases  is  determined  by  the  compo- 
sition,  temperature 5  and  pressure  of  the  medium,  and  the  wave  length 
of  the  radiation. 

The  use  of  the  exponential  attenuation  function  generally  requires 
the  solution  of  very  complicated  integral  relationships  to  obtain  the 
radiative  properties  of  a  given  medium.  A  volume  of  gas  participating 
in  a  radiant  interchange  emits  from  its  surrounding  envelope  or  absorbs 
within  itself  an  amount  of  energy  equal  to  the  sum  of  the  contributions 
from  each  differential  volume  element.  Functions  describing  the  location 
of  differential  elements  within  an  absorbing  region  are  generally  complex 
and,  when  incorporated  in  radiation  equations,  are  not  readily  integrated. 
One  simple  case,  however,  is  that  of  a  small  surface  at  the  center  of  the 
base  of  a  gaseous  hemisphere.  The  exchange  of  radiation  between  the 
surrounding  gas  and  the  surface  element  is  readily  calculated,  since  the 
distance  to  the  hemispherical  surface  is  constant  and  is  independent  of 
direction. 

12  1 

Hottel  and  co-workers  ^  ^  noted  that  a  given  system  had  an 
"equivalent”  dimension  (or  mean  beam  length)  such  that  a  hemisphere  of 
gas  having  this  "equivalent"  dimension  for  its  diameter  would  radiate  to 
a  differential  area  at  the  center  of  its  base  at  a  rate  equal  to  that  of 
a  small  area  on  the  surrounding  surface  of  the  real  gas  shape.  Then,  if 
the  "equivalent"  dimensions  were  known,  practical  problems  involving  the 
characteristics  of  the  radiation  could  be  calculated  from  the  simplified 
case.  Hottel  selected  the  most  frequently  occurring  furnace  shapes  and 
calculated  the  "equivalent"  dimensions.  The  calculations,  unfortunately, 
had  to  be  further  restricted  to  arbitrarily  selected  dimensions  because 
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the  results  were  dependent  on  the  radiation  attenuation  coefficients 
of  the  absorbing  medium.  Results  were  obtained  for  certain  cases  where 
the  coefficient  could  be  approximated  as  being  zero  and  for  other 
instances  where  finite  coefficients  are  frequently  encountered  in 
practice;  other  values  not  specifically  evaluated  were  recommended  as 
being  around  85  per  cent  of  the  result  for  the  zero  coefficient. 
Although  limited,  these  calculations  sufficed  until  recently  for  the 
majority  of  the  radiation  calculations  involving  absorbing  gases. 

The  next  step  was  to  extend  the  method  to  apply  to  particle 
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clouds.  Wollenberg  and  Morrow  and  Haslam  and  Hottel  considered  the 
case  where  energy  was  removed  from  a  radiant  beam  by  the  screening 
effect  of  randomly  oriented  particles.  Based  on  probability  consider¬ 
ations,  the  relationship 
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was  derived  which,  except  for  the  absorption  coefficient  k,  was  iden¬ 
tical  to  the  Beer's  law  function  for  gases.  The  constant  for  the 
modified  equation  was  the  sum  of  the  particle  cross-sectional  areas 
divided  by  the  sum  of  the  particle  volumes  times  the  particle  volume 
concentration;  for  monodispersed  systems  and  for  dispersions  considered 
to  be  represented  by  an  average  particle  size,  the  constant  reduces  to 
the  average  particle  cross-section  times  the  particle  number  concen¬ 
tration.  This  exponential  relationship  was  applied  in  calculating  the 
emissivity  or  absorptiArity  of  a  pulverized  coal  flame  after  assumptions 
were  made  about  the  particle  size  distribution  and  the  rate  of  combus¬ 
tion  of  the  individual  particles .  It  was  noted,  however ,  that  the 
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results  were  almost  inevitably  lower  than  measured  flame  emissivities. 
The  discrepancy  was  attributed  to  the  contribution  of  cracked  hydro¬ 
carbons  producing  luminous  as  well  as  residual  ash  particles  not 
considered  in  the  theoretical  derivation.  The  problem  was  not  pursued 
further  because  most  pulverized  coal  installations  involved  such  large 
flames  that  the  emissivities  did  not  differ  greatly  from  unity. ^ 

Other  calculation  methods  were  devised  by  Johnstone,  Pigford, 
and  Chapin.”^  These  dealt  with  the  heat  transfer  to  a  cloud  of  parti¬ 
cles  falling  through  a  heated  pipe.  The  objective  was  to  provide 
information  for  use  in  the  "flash"  heating  of  finely  divided  particles 
by  passing  them  through  a  section  of  pipe  at  temperatures  as  high  as 
1050°  F.  The  carbonization  of  pulverized  coal,  the  briquetting  of 
blast  furnace  dusts,  and  the  production  of  soluble  phosphates  afforded 
practical  applications  for  the  results.  Convection  was  the  dominant 
mode  of  heat  transfer.  For  wall  temperatures  near  1000°  F,  however, 
the  particle  clouds  gained  up  to  approximately  20  per  cent  of  their 
energy  by  radiation.  The  small  absorption  of  radiation  was  due  to  the 
large  size  of  the  particles  employed,  which  were  sand  grains  on  the 
order  of  500  to  1000  microns  in  diameter.  To  calculate  the  radiant 
heat  transfer,  the  investigators  took  advantage  of  Hottel's  basic 
attenuation  relationship  but  used  developments  by  Nusselt  and  Jakob 
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to  evaluate  particle  cloud  emissivities.  Combined  analytical  and 
graphical  procedures  were  employed  to  integrate  modified  relationships 
and  to  obtain  results  for  a  sphere,  an  infinite  cylinder,  and  an  in¬ 
finite  parallel  plate.  Generalized  curves  were  then  constructed  as 
shown  in  Figure  1,  and  approximate  radiation  properties  for  particle 
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Particle  Cloud  Absorption  Efficiency, 


2kR 

Figure  1.  Absorptivity  of  Clouds  of  Particles  in  Furnaces 
of  Various  Shapes. 


clouds  were  estimated  from  qualitative  deductions  about  the  relative 
absorptance  of  the  three  different  cloud  configurations.  For  a  given 
value  of  2kRj  a  particle  cloud  between  infinite  parallel  plates  was 
noted  to  have  the  greatest  absorptivity.  Qualitatively  this  was 
expected  since  the  average  length  of  a  radiant  beam  should  be  largest 
for  infinite  parallel  plates  because  of  the  two  infinite  directions. 

The  average  beam  length  should  be  intermediate  for  a  cylinder  with 
one  infinite  direction  and  least  for  the  sphere  with  no  unlimited 
dimensions . 

The  investigation  concluded  that  the  absorptivities  of  cloud  con¬ 
figurations  encountered  in  practice  could  be  estimated  by  interpolation 
between  the  curves  of  the  idealized  cases.  For  instance^  a  cubic  cloud 
structure  should  have  an  emissivity  between  that  of  a  cylinder  and  a 
sphere  while  a  long  rectangular  enclosure  should  result  in  a  value 
between  the  cylinder  and  parallel  plates.  The  method  thus  afforded  a 
means  for  estimating  particle  cloud  emissivities;  it  was  not  intended 

to  give  precisely  an  evaluation  of  a  cloud’s  shape, 
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In  a  later  investigation,  Oktay^  extended  this  method  of  attack 
to  include  particles  as  small  as  50  microns.  As  would  be  expected,  the 
radiation  contribution  to  the  overall  heat  transfer  increased  rapidly 
with  decreasing  particle  size.  Shulderberg,  et  made  a  similar 

investigation  of  the  heat  transfer  characteristics  of  solid-gas  suspen¬ 
sions  for  possible  use  as  the  working  fluid  of  a  nuclear  power • genera¬ 
tion  system.  Again,  high  rates  of  heat  transfer  were  obtained  for  very 
small  particles.  This  latter  investigation  was  exploratory  in  nature 
and  did  not  emphasize  basic  analyses. 


Hoffman  and  Gauvin^^  conducted  a  more  fundamental  study  in  which 

the  evaporation  of  a  spray  contained  in  a  high- temperature,  cylindrical 

column  was  investigated.  This  analysis  of  radiant  heat  transfer  to  a 
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cloud  of  particles  amounted  to  an  extension  of  the  Hottel-Cohen  method 
in  which  the  system  was  divided  into  zones,  and  the  radiation  between 
these  zones,  expressed  as  geometric  interchange  factors',  was  evaluated. 
Results  for  the  total  system  were  obtained  from  energy  balances  over 
each  zone.  The  zones  consisted  of  portions  of  the  enclosing  surface  and 
of  the  absorbing  medium  that  were  small  enough  in  size  to  be  considered 
isothermal.  To  be  exact,  the  zones  must  be  of  differential  size,  yet, 
because  of  the  lengthy  computations  involved,  the  number  of  zones  were 
restricted.  Hoffman  and  Gauvin  compromised  by  taking  cylindrical  aero¬ 
sol  volimes  of  height  equal  to  diameter  and  by  developing  calculations 

for  three  such  volumes.  Mean  beam  lengths  as  previously  outlined  by 
11 

Haslam  and  Hottel  are  derivable  from  this  work. 
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In  a  more  recent  analytical  study,  Einstein  continued  the 
zone-type  of  calculation  and  made  a  computer  analysis  of  the  heat 
transfer  to  an  absorbing  medium  within  a  black  circular  pipe  under  the 
combined  influence  of  radiation,  gas  flow  through  the  pipe,  and  thermal 
conduction  from  the  wall  to  the  gas.  Here,  the  zones  consisted  of  short 
lengths  of  cylindrical,  concentric  shells,  and,  at  the  outer  boundary, 
of  the  enclosing  portion  of  the  emitting  pipe  surface.  Specific  results 
were  obtained  for  a  pipe  with  absorbing  medium  opacities  of  1.0  to 
10.0  [opacity  being  the  absorption  coefficient  times  the  diameter  of 
the  pipe  2R]. 


9 


15 

Lanzo  and  Ragsdale  performed  an  experimental  study  of  radiant 
heat  transfer  to  a  particle  cloud  flowing  in  a  quartz  annulus  surround¬ 
ing  an  electric  arc.  Very  small,  submicron  (altho\:igh  only  partially 
dispersed)  were  shown  to  attenuate  considerable  quantities  of  energy 
and  to  make  excellent  radiation  absorbers.  Even  better  results  would 
undoubtedly  have  been  obtained  had  good  particle  deagglomeration  been 
obtained. 

All  the  foregoing  studies  were  similar  in  that  the  basic  objective 
was  to  evaluate  the  net  heat  interchange  between  a  solid  surface  and  a 
particle  cloud  provided  the  absorption  characteristics  of  the  cloud  were 
known.  Unfortunately,  many  liquid  sprays  and  particle  dispersions  may 
be  highly  reflective  or,  because  of  particle  shape  and  dimensions  and 
size  distribution,  exhibit  very  complicated  radiation  scattering  pheno¬ 
mena  and  as  a  result  possess  uncertain  radiation  properties. 

An  entirely  different  approach  is  required  to  account  for  these 
parameters.  Early  work  of  the  latter  type  attempting  to  define  the 
emissive  properties  of  luminous  flames  in  furnaces  is  described  in 
Heat  Transmission  by  McAdams. Semi-empirical  results  are  represented 
by  the  equation 
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where  is  the  spectral  emissivity  of  a  particle  cloud  of  concentration 

C  with  a  mean  beam  length  L  and  the  function  f(l)  is  an  absorption 
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coefficient.  The  function  is  represented  either  by  (k/A  ‘  ,  -which  is 

supposed  to  account  for  particle  size  variations,  or  by  [l  -  15(r/x)  ]/X, 
where  r  is  the  radius  of  a  specific  particle.  The  general  validity  of 
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these  functions  is  uncertain however.  Stull  and  Plass  attacked  the 

problem  more  fundamentally  by  using  spectroscopic  and  electromagnetic 
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data  to  solve  the  Mie  equations  for  the  absorption  and  reflective 
characteristics  of  small  luminous  particles.  Based  on  observed  parti¬ 
cle  size  distribution  data,  emissive  properties  of  luminous  carbon 
flames  were  obtained  for  a  temperature  of  2250^  K  and  exposure  to 
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several  appropriately  selected  radiation  wave  lengths.  Thring,  ^  al. 
applied  the  data  of  Stull  and  Plass  to  practical  cases  and  found  that 
for  flames  two  meters  across  (square  tunnel)  the  emissivity  could  be 
predicted  within  ±  15  per  cent,  but  for  smaller  flames  of  0.5  meter 
across  deviations  of  about  ±  300  per  cent  were  observed.  Strangely, 
small  laminar  flames  in  the  laboratory  had  emissivities  in  good  agree¬ 
ment  with  the  theoretical  predictions. 

Ordinarily,  solid  particles  must  be  on  the  order  of  one  micron 

or  less  to  exhibit  pronounced  scattering,  almina,  quartz,  and  other 
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visually  transparent  minerals  being  exceptions.  ^  Liquids,  on  the 
other  hand,  have  significant  transmissive  and  refractive  properties. 
Because  water  and  atomized  fuel  droplets  occur  frequently  in  high  tem¬ 
perature  processes,  investigations  of  their  radiation  properties  have 
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been  made.  Thomas  determined  the  absorptivity  of  water  droplets  by 
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use  of  geometric  optics.  Friedman  and  Churchill  produced  similar 

absorptivity  results  for  fuel  droplets.  The  data  in  either  case  are 

only  valid  for  droplets  larger  than  the  radiation *s  wave  length.  The 

error  is  uncertain  if  this  condition  is  not  met  since  the  optical  laws 

of  reflection  for  large  spheres  do  not  hold  for  bodies  smaller  than  the 
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radiation  wave  length. 
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still  another  problem  associated  -with  the  radiant  heating  of 
particle  dispersions  is  the  dissipation  of  the  energy  absorbed  by  the 
particles.  Mathematical  solutions  have  been  obtained  for  the  transfer 
of  heat  from  submerged  bodies  to  a  surrounding  medium,  but,  because 
of  the  many  complicating  factors,  they  were  highly  idealized.  Heat 
transfer  in  the  neighborhood  of  a  particle  is  dependent  upon  the 
particle  shape,  size,  distance  to  other  particles,  and  the  time-depen- 
dent  hydrodynamic  flow  patterns  in  the  surrounding  space.  Obtaining 
accurate  descriptions  of  all  the  variables  required  for  an  exact 
analysis  has  not  been  feasible.  With  certain  simplifications,  however, 
guidelines  have  been  established  for  estimating  the  heat  transfer 

between  suspended  particles  and  their  surroundings. 
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Notably,  Sliecher  and  Churchill  investigated  the  problem  for 
quiescent  aerosols  and  derived  limiting  results  where . spherical,  uni¬ 
formly  dispersed  particles  are  heated  by  isotropic  radiation  and  cooled 
by  conduction.  The  development  was  made  first  for  a  single  particle  in 
an  infinite  medium  and  was  then  extended  to  apply  to  a  dense  dispersion 
of  particles.  For  a  single  particle  in  a  large  medium,  its  temperature 


history  was  found  to  be  expressed  exactly  by 


t  -t 
0  o 


i  ■] 


c»  T  ^  2  /  ^  2 

2  r  /PpOp»  ^  1 

;;  J  , 

l-(2f-l)z  +$  z 


(1.4) 
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the  latter  is  not.  At  high  values  of  §,  however,  the  two  become  nearly 
identical.  Since  this  number  is  usually  large  (over  Uoo),  both  equa¬ 
tions  in  practical  applications  were  found  to  agree  within  ±  3  per  cent. 
The  functions  were  plotted  in  non-dimensional  form  as  shown  in  Figure  2 
for  em  arbitrary  vailue  of  $  =  500,  which  corresponded  to  a  particular 
case,  viz.,  a  25  micron  coal  particle  exposed  to  a  2000°  F  black-body 
heat  source.  For  these  conditions  a  dimensionless  abscissa  of  4  repre¬ 
sents  an  exposTire  time  of  l6.1  milliseconds  during  which  the  coal 
particle  attains  over  95  per  cent  of  its  maximum  possible  temperature. 

For  a  dense  particle  dispersion,  the  ten^erature  rise  of  a  single 
particle  could  be  described  by  the  function 


t 


-t 
0  ’^o 


+ 


4TTank^9 

~lc 


■) 


(1.6) 


containing  two  linear  terms  and  one  e3q)onential  term.  Initially  the 
exponential  term  is  controlling  and  a  particle’s  temperature  rises  in  a 
manner  very  nearly  that  of  a  single  particle  in  an  infinite  medium. 

This  characteristic  is  demonstrated  graphically  in  Figure  3>  where, 
below  an  abscissa  of  1.0,  the  curves  shown  are  very  nearly  exponential. 
For  the  coal  particle,  an  abscissa  of  1.0  again  represents  l6.1  milli¬ 
seconds.  At  still  greater  times,  the  particle  temperature  establishes 
a  pseudo- equilibrium  condition  such  that  the  energy  received  by  radia¬ 
tion  equals  the  heat  lost  by  conduction;  the  linear  terms  are  then  con¬ 
trolling  as  may  be  noted  from  its  long  term  temperature  history  as 
shown  in  Figure  4. 

The  temperature  difference  at  pseudo- equilibrium  between  the 
particles  and  the  surrounding  air  is  constant,  and  this  difference  for 


(t_-t  )/aqa/k.  =  Non-dimensi 

/aqa/k^  =  Non-dimensi bnali zed  Temperature  Ordinate 

Temperature  Ordinate 


Radiation 


the  example  is  3l8^  F.  Between  the  initial  and  pseudo-equilihriiim  con¬ 
ditions  there  is  an  intermediate  range  where  all  the  terms  are  signifi¬ 
cant.  This  behavior  is  controlled,  as  shown  in  Figures  3  and  4,  by  the 
parameter,  4n^^C^p^C^/3pG5  which  represents  the  ratio  of  the  particle 
heat  capacity  to  that  of  the  air  associated  with  each  particle. 

The  calculations  of  Sleicher  and  Churchill  are  important  because 
they  indicate  limiting  conditions  and  temperature  differences  that  may 
be  expected  in  an  actual  system.  The  significant  finding  is  that  time 
intervals  of  the  order  of  milliseconds,  and  sometimes  micro-seconds, 
are  required  for  essentially  corrp)lete  transfer  to  occur.  In  flowing 
systems  this  means  th  t  only  very  short  distances  are  required  for 
particle-gas  interactions.  For  example,  a  gas-solid  suspension  flow¬ 
ing  at  50  feet  per  second  through  a  pipe  and  receiving  2000^  F  black 
body  radiation  would,  if  the  particles'  were  50  microns  or  less,  estab¬ 
lish  pseudo-equilibri-um  after  only  6  inches  of  movement  along  the  pipe, 
according  to  the  predictions.  After  leaving  the  heated  pipe,  the  sus¬ 
pended  particles,  would  in  even  less  distance  return  very  nearly  to  the 
ambient  gas  temperature. 

Although  the  above  type  of  thermal  reaction  is  rapid,  even  less 
time  than  predicted  is  frequently  required.  Since  the  development  was 
for  quiescent  aerosols,  the  particles  were  considered  stationary  rela>- 
tive  to  the  gas,  and  conduction  only  was  assumed  to  exist  throughout  the 
region  among  the  particles.  In  flowing  systems,  however,  there  must 
be  convective  transfer  and  conduction  can  only  exist  in  a  small  boundary 
layer  surrounding  a  particle.  Conduction  due  to  the  decreased  thick¬ 
ness  is  actually  faster  and  the  temperature  of  a  particle  more  quickly 
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adjusts  to  that  of  the  surrounding  air.  Under  highly  accelerated 
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particle-gas  conditions,  however,  slippage,  i.e.,  the  particle 
lagging  in  velocity  behind  that  of  the  surrounding  air  because  of  momentum 
differences,  can  change  a  particle's  environment.  This  would  produce  the 
opposite  effect  on  the  difference  between  particles  and  average  air  tern- 
perature.  Carlson  measured  spectroscopically  the  average  air  and 
particle  temperatures  in  a  rocket  engine  where  severe  acceleration  was 
exhibited.  There,  magnesium  oxide  particles  on  the  order  of  one  micron 
lagged  the  gas  temperature  about  650^  F  when  the  rocket  exhaust  was  about 
4200^  F.  This  temperature  difference  is  considerably  larger  than  would 
be  expected  from  Sleicher  and  Churchill's  calculations.  Where  large 
accelerations  are  e^^pected,  the  effect  must  be  investigated  further, 
therefore. 

That  heat  transfer  from  a  particle  to  its  surroundings  is  rapid 
may  be  readily  demonstrated  by  considering  the  classical  developments 

28 

for  heat  transfer  from  a  sphere  submerged  in  a  gas.  The  limiting  case 
is  for  a  Wusselt  number  of  2  where  only  conduction  prevails.  Thus,  the 
lowest  heat  transfer  coefficient  is 


h 


a 


(1.7) 


If  the  macroscopic  properties  of  the  gaseous  medium  are  not  violated,  a 
small  diameter  particle  has  associated  with  it  a  very  high  heat  transfer 
coefficient,  even  for  conduction.  For  a  100  micron  particle  in  air,  h^ 
is  over  9000  Btu/hr  ft^  °F  as  normally  encountered  in  convective  heat 
transfer.  As  may  be  seen  from  equation  (l.?),  the  increased  heat  tranfer 
coefficient  arises  from  its  inverse  variation  with  particle  diameter.  Use 


17 


has  been  made  of  this  fact  in  measuring  air  temperatures  in  the  presence 
29 

of  radiant  fields.  When  the  air  temperature  is  determined  with  several 
thermocouples  of  successively  smaller  size^  the  actual  and  indicated  tem¬ 
peratures  become  more  nearly  equal.  The  results  are  extrapolated  to  zero 
diameter  to  obtain  the  actual  gas  temperature. 
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CHAPTER  II 


THEORETICAL  INVESTIGATIONS 
Method  of  Analysis 

It  is  frequently  expedient  to  approach  radiation  problems  by 
first  assuming  black  body  behavior  for  all  the  bodies  involved,  then, 
if  a  solution  is  obtained,  to  seek  progressively  more  accurate  des¬ 
criptions  of  the  problem  by  modifying  the  basic  result  with  the  para¬ 
meters  necessary  to  account  for  the  non-ideal  characteristics  of  the 
system.  The  fundamental  developments  of  classical  radiative  heat 
transfer  theory  have  followed  this  procedure.  Accordingly,  the  heat 
exchange  between  a  particle  cloud  and  a  cylindrical  enclosure  was 
first  investigated  when  both  were  perfect  radiators.  A  combined 
analytical  and  numerical  approach  was  employed  to  derive  an  accurate 
evaluation  of  the  radiant  heat  transfer  to  a  particle  cloud  within  an 
enclosure  having  practically  any  length-to-diameter  ratio.  Calculated 
results  are  readily  obtained  by  use  of  a  digital  computer. 

Development  of  Equations 

Radiation  emitted  from  each  of  two  black  surfaces  of  differential 
size  and  received  by  the  other  is  described  by 


■^<12  =  Snl  h  *"2  “1 

(2.1) 

and 

<21  ■  ‘bn2  ®2 

(2.2) 
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for  which 


(2.3) 


2 


cos  P 

-:x-^2 


cos  3 

*■1  '  -^2-  '“l  ("•'*) 

The  subscript  b  indicates  black  body  behavior  and  n  specifies  that 
the  radiation  intensity^  i,  is  normal  to  the  emitting  surface.  The 
number  subscripts^  1  and  2,  refer  to  the  surfaces  exchanging  radiation. 
Figure  5  defines  the  geometry  of  the  plane  angles,  3^^  Q'^d  3^3  the 
solid  angles,  dou^  and  dco^. 

Since  the  surfaces  are  black,  all  the  impinging  radiation  is 
absorbed.  The  net  exchange  of  energy,  therefore,  is 


‘^S321  "  '^‘^12 

If  equations  (2.3)  and  (2.4)  are  substituted  into  (2.l)  and  (2.2)  and 
the  results  applied  to  (2.5)?  the  net  heat  transfer  is 

cos  B.  cos  8  V 

^\n2"^bnl^  (  72  /  (2.6) 

b 

or  in  integrated  form 

P  r*  .cos  R  cos  &  K 

<-J  J  (2.7) 

A, 

Equation  (2.7)  is  a  general  expression  that  relates  the  radiant  heat 
transfer  between  two  finite  surfaces.  For  a  specific  application,  the 
spatial  arrangement  of  dA^  and  dA^  in  the  differential  expression  must 
be  described  with  appropriate  functions  to  permit  the  integration  indi¬ 
cated  in  equation  (2.7).  For  a  particle  cloud  receiving  radiation  from 
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(a)  Orientation  and  Separation  of  Differential  Areas 


(b)  Solid  Angle  Relationships 


Figure  5.  Angular  Functions  Used  in  Defining  Radiation  Interchange 
Between  Two  Areas  of  Differential  Size. 


a  tubular  enclosure ^  adaptation  was  made  as  indicated  in  Figure  6  where 
the  system  of  coordinates  is  also  defined.  As  shown,  radiation  leaves 
the  tube  wall  with  an  intensity,  i^,  and  proceeds  along  the  line  segment, 
S,  to  a  small  elemental  volume  of  aerosol.  The  location  of  dA^  is  des¬ 
cribed  by  the  variables  y,  z  and  Y,  and  the  position  of  the  elemental 
volume,  by  z  and  r.  Other  quantities  of  equation  (2.7)  in  terms  of 
these  coordinate  variables  are: 


dA  =  RdH'dy 
w 

(2.8) 

o  722 

S  =  Vm  +y 

(2.9) 

/R^+rl(2Rr)cos  Y 

(2.10) 

and 


cosp 

w 


S^+R 


2 

-r 

2RS 


2 


2 


(2.11) 


Receiving  the  radiation  is  dA  which  is  the  surface  area  of  the 

s 

particles  contained  in  the  volume  element.  As  shown  in  Figure  how¬ 
ever,  the  effectiveness  of  an  area  fer  intercepting  radiation  is  deter¬ 
mined  hy  its  projection  perpendicular  to  the  incoming  radiation.  This 
projection  is  obtained  from  multiplication  of  the  area  by  cos  3^  which 
is  the  cosine  of  the  angle  between  the  normal  direction  to  the  surface 
and  the  impinging  radiation.  When  a  particle  is  considered,  the  pro¬ 
jection  is  its  cross-sectional  area,  A^.  The  effective  particle  area 
participating  in  a  radiant  interchange,  then  is  given  by 

dA  =  CS  dV  (2.12) 

c^  V 
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End  View 
Projection 

Area  on  Wall 

Longitudinal  Pictoral  View 

6.  Definitions  of  Geometric  Variables  for  Tube  Radiating 
to  Enclosed  Particle  Cloud. 


or 

dA  =  CS  rdrdGdz  (2.13) 

c  V 

where  is  the  cross-sectional  area  per  unit  volume  of  powdered  material 
(ft^/ft^),  and  C  is  the  volume  concentration  of  powder  in  the  aerosol 
(ft^  solids/ft^  aerosol).  The  quantity  may  be  computed  from  statis¬ 
tical  measurements  of  randomly  oriented  particles.  If  the  surface  area 
of  a  powder  is  known  from  other  information,  the  required  area  is  0.25 
times  the  envelope  surface,  or 


dA  =  0.25  dA 
c  s 


(2.14) 


Equation  (2.l4)  follows  from  the  fact  that  the  projected  area  of  randomly 

oriented  particles  whose  surfaces  are  everywhere  positively  curved  is 
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0.25  times  the  surface  area. 

Substitution  of  the  required  quantities  into  the  radiation  equa¬ 


tion  yields 


2n  R  2tt  L  L-z 

f  f  (i. 

J  J  bnw  bnp 
o  o  o  o  -z 


)  (CS^)  (i)  )  ( r )  dydzdf  drde  (2.15) 

s 


where  the  limits  of  integration  are  as  shown. 

The  result,  however,  is  still  incomplete.  Radiation  passing 
through  a  cloud  has  its  intensity  attenuated.  For  black  particles  this 
reduction  may  be  e^qpressed  as 

i  =  i  e"^^  (2.16) 

s  w 
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where  k,  is  a  function  of  particle  geometry  and  concentration;  the 

radiation  intensity  i  at  the  wall  is  reduced  to  i  after 

w  s 
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Solution  for  Integrals  of  Heat  Transfer  Equation 
The  indicated  integrations  remain  to  be  performed  after  all  the 
quantities  of  equation  (2.17)  are  evaluated.  The  integral  expression 
may  be  simplified  by  considering  the  particles  not  to  radiate  back  to 
the  wall  either  by  radiation  scattering  or  by  significant  thermal  radia¬ 
tion.  The  development  is  then  valid  for  the  radiant  heating  of  aerosols 
by  direct  transmission,  which  is  of  prime  interest  in  this  investiga¬ 
tion.  Secondary  emissions  and  reverse  thermal  radiation  is  more 
appropriately  considered  in  a  separate  analysis  since  cannot  be 

taken  as  constant  and  be  removed  from  the  integration  process.  It  must 
be  further  assumed  that  the  radiation  intensity  of  the  emitting  tube  wall 
is  constant  over  the  entire  surface.  These  modifications  yield  for  the 

total  heat  absorbed 

2tt  R  2n  L  L-z  -kS  2  2  2  2 

II  I  (2.18) 

o  o  o  o  z 

A  second  change  may  now  be  made  to  treat  the  integration  with  respect  to 
the  first  two  variables,  y  and  z,  as  one  operation  which  is  a  lengthy 
development  itself  and  then  to  integrate  the  remaining  ones  separately. 
The  following  notation  will  be  employed 


2tt  R  2n 


R  2tt 


where 


L  L-z 


I  =  J  J  (^)( 


-k  „2  „2  2  2 

^/S  +R  -r  -y 


)dydz  . 


(2.20) 


The  9  integration  is  performed  directly. 

The  evaluation  of  I  is  initiated  with  an  expansion  of  the  expo¬ 


nential  term.  Thus, 


2  2  ^  4  4 

e'^  =  1  -  kS  +  —  -  ^  ^  - 


(2.21) 


The  cos  B  factor  can  also  be  expressed  as 
w 

S^+R^»r^-y^  =  ^  +  Sfzl?- 

2RS  2RS  2RS  2RS 


(2.22) 


With  these  changes  equation  (2.20)  becomes 


L  L-z 


2  R  2  4  R 


o  -z 

L  L-z  ^  2  22  “^2  42 

-  (a)  J  J  -  5-  *  M-  -  ^ 

^  S  ^ 

o  -z 

2  2  1/2 

In  summation  notation  and  with  (m  +y  )  '  '  substituted  for  S  the  terms 
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I  =  I  +  (2.24) 

O'  3 

where  and  3^  are  integrals  of  the  two  suiranations. 

Integrations  for 

Because  of  the  mathematics  involved,  it  is  necessary  to  perform 
the  integration  in  four  separate  divisions  classified  according  to  the 
position  of  the  term  in  the  series  being  integrated.  Two  classifications 
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are  simply  for  the  first  and  second  terms  of  the  series,  a  third  for 

all  successive  odd  terms  and  a  fourth  for  all  subsequent  even  terms.  The 
terms  are  even  or  odd  according  to  the  value  of  n.  The  total  integral, 
therefore,  is 


I 

O' 


+  I  +I-  +  I) 
cvl  Q'2  aS  0^4 


(2.25) 


where 

0^3  O' (odd) 


^^(even) 


Evaluation  of  and 


The  initial  term  of  the  series  is 


L  L-z 


J  J  [( 2/2,372] 

o  -z  ^  y 


(2.26a) 


which  becomes 


= (3)  j  [■ 


m  [m^+(L-z)^]^/^ 


]  + 
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[7™]} 

[m  +z  ]  ’ 


dz 


and 


2  2 


^.1 = <5)  - "] 
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(2.26b) 


(2.26c) 


The  second  term  of  the  exponential  expansion  is 
L  L-z 


^02  -  (S)  J  J  (■>‘)('>"-'^^)[72S7] 

(m  +y  ) 


0  -z 


which  by  successive  integration  and  simplification  becomes 
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and 


I  =  (-k)(R  -r„)  r  -1  (Lr£)  +  tan"^  (-)]  dz  (2.27b) 
a2  2mR  J  ^  m  \  i  / 

o 

^02  ‘  [(2L)tan-l  (i)  -  (m)  In  (1  +  4)]  (2.27c) 

In  the  integration  of  equation  2.27b  use  is  made  of  the  trigonometric 
identity^ 

tan  ^  (— )  =  -  tan  ^  (— )  (2.28) 


Also 5  since  z  and  m  are  both  positive  and  their  ratio ,  z/m^  varies  from 
zero  to  infinity,  the  inverse  tangent  of  (L/m)  involved  only  the  prin¬ 
ciple  values  of  the  function. 

Evaluation  of  I  Since  n  is  odd  and  larger  than  1,  the  exponent 
_ Q!3 

(— ^)  in  equation  (2.23c)  is  a  member  of  the  set  of  fractions,  -l/2,  l/2, 
3/2,  5/2,  7/2,  9/2,  etc.;  and  the  complete  integration  for  the  first  n 


terms  is 


L  L-z 


=  <S)  I  {  ^  }  j  J  (2.2, a) 


n=3 

odd 


o  -z 


For  a  particular  value  of  n.  the  relationship  reduces  to 

^c3n  =  (a>  { [i5^]  *  [isttF-''^]  j  / 


o  -z 


This  equation  then  can  be  integrated  by  the  following  expressions  from  a 
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table  of  integrals: 


I 


2cx+b  n  7 —  2n+l 

4(n-i-l)c  ^  ^  2(n+l)f 


dx  (2.31) 
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J 


dx 


/x 


'/x  dx 


(2.32) 


and 


I 


1  ^ 

Vl(  '/"c 


sinh"^  (_f£^L_)  if  c  >  o  (2.33) 


4ac+b' 


where  x  =  +  hx  +  a, 

4c 


and  f  = 


2  • 


4ac-b 

Equation  (2.3l)  is  more  conveniently  employed  when  transformed  to  yield 


r  p  .  2cx+h  p  2(p-'5')+1 

=  "2Mr 


^ -  dx 


(2.34) 


The  integral  formulas  can  he  applied  to  equation  (3.30a)  hy  making 

2  2  2 
X  =  m  +y  .  Then  it  follows  that  x=y,  a=m,b=0,  c=l,  and 

f  =  l/m^. 

The  first  two  members  of  series  are  integratable  directly.  The 
results  are 

^.33  -  {  [w]  ^  }  1  <i)] 

O 

2  ^ 

=  (^)  1  +  (^)(R^-3^^)J  j  [  Sinh"^  (^^)  +  Sinh  ^  (|)J  dz  (2.30b) 


and 
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^.33  =  <a)  { [« *  [TfiTTiF-^ }  J  { [¥F^(-^)^T  ^ 

O 

|][m^+z^]^  +  (^)  +  sinh"'^  (|)]  j  dz  (2.30c) 


dz 
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(a^t 


“  ‘‘-  +  (t|-){1>^-i'^)]  J  [si*"'’  (^)  +  sliih'^  (|)]'^^  (2-30c) 


where  the  final  subscripts  3  and  5  refer  to  the  particular  value  of 

n. 

By  use  of  equations  (2.30a)  and  2.34)  the  integral  becomes  for 
higher  values  of  n 

L  T 
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y  M,Y 


} 


(2.30d) 


where  p  is  a  member  of  the  set  of  exponents  and  the  notation  (MjY)^ 

2  2 

symbolizes  the  integrand,  (m  +y  )  2  .  A  second  application  of 

equation  (2.3^)  gives 
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(2.30e) 
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and  repeated  use  of  the  integration  formulas  produces  the  series 


51 


=  (s*  1  [iferr] 
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(M,Y)' 


(2.30f) 


The  series  continues  until  the  exponent  (M,Y)  in  the  last  integral  is 
(-I/2),  and  for  this  case  equation  (2.33)  is  employed  to  integrate  the 
final  term.  With  the  substitution  of  for  p,  the  complete  result 

for  any  n  value  greater  than  5  can  he  written 

/I  V  r  r  1  .  r  t„2  21 1  r  J/  y  2._2x  ^ 
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and  with  the  limits  evaluated,  it  becomes 
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(2.30h) 


52 


Next  to  t»e  performed  is  the  integration  with  respect  to  z,  which 

Q  2  2  2 

involves  power  functions  of  [m  +  (L-z)  ]  and  [m  +  z  ]  and  the  inverse 
hyperbolic  relations,  tinh"^  {—)  and  sinh"^  (|) .  The  former  are 
treated  directly  hy  making  the  substitutions ^ 

2  fr 

u  =  m  +  (Ij-zj 

du  =  -2(L-z)dz 

2  2 
V  =  m  +  z 

dv  =  2zdz, 

and  integrating.  The  result  is  still  a  set  of  similar  power  functions. 
The  inverse  hyperbolic  sines  are  collected  and,  when  integrated  by  parts, 
become 

L  I 

J  +  sinh'^(|)]dz  =  2m  +  (2L)  sinh‘^(i)  .2im^+L^)  (2.34) 

o 

Equations  (2.30b),  (2.30c)  and  (2.30h),,  after  integration  with  respect  to 
z  are  then  summed  for  all  values  of  n  according  to  equation  (2.29a) 
to  yield  the  final  result 
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j^(2m)  +  (2L)  sinh“^(^)  -  2(L^+m^)  J  j- 

Evaluation  of  The  expression  for  I  ,,  is 
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Since  n  is  even  and  larger  than  2,  the  exponent  (n-4)/2  is  an  integer. 
The  quantity  to  he  integrated  then  could  be  expressed  with  the  binomial 
theorem  as 

(F) 
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Likewise  the  integration  with  respect  to  z  can  be  readily  accomplished, 
the  final  result  being 
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w  } 

Integrations  for  (See  equation  2.24) 

The  evaluation  of  the  function  is  very  similar  to  that  for 


The 


result  is  somewhat  different,  however,  because  of  the  y  term  in 


the  integral  which  is 

L  L-zn=n 


I  =  -(— ) 
^2R'' 


o  -z  n=l 


As  before,  it  is  necessary  to  treat  the  first  two  terms  individually. 

The  remaining  ones  can  then  be  integrated  by  parts.  This  operation  pro- 

duces  a  set  of  integrals  exactly  like  those  for  ^Qf4 

the  coefficients.  The  integrations  for  these  situations  are  therefore 

combined. 
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Evaluation  of  and  I„_.  The  expression  for  I.,  is 
_ 31 _ B2  ^  pi 


L  L-z  2 

J  J  /  2  ^2,.3/2 

o  -z  ^  ' 


(2.38a) 


which  by  use  of  a  formula  from  a  standard  set  of  integral  tables  becomes 
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+  sinh 


(2.38b) 


The  terms  containing  square  roots  integrate  directly,  while  the  parts 
formula  was  required  for  the  inverse  hyperbolic  sines.  After  inter¬ 
mediate  operations  and  simplifications,  the  final  result  is 

1/2 

^B1  +  (L)  sinh"^(^)J  (2.38c) 


The  second  term  of  the  series,  Ipp?  is 

L  L-z 


=  -  J  J 

„  „  (m  +y  ) 


(2.39a) 


The  evaluation  of  this  term  is  straightforward  and  almost  like  the  one 
above.  The  results  for  the  two  integrals  are 


^32 


J  +  z  -  (m)'fcan"^(4)l  dz  (2.39^) 
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^B2  [L^-(2inL)tan"^(^)  +  (m^)  ln(l4^)  J 


(2.39c) 


I  Terms  for  n  >  2.  The  general  expression  for  I  is  given  hy 
_  P 

equation  (3*37a),  With  n  larger  than  2  the  relationship  can  be 
integrated  by  parts,  the  parts  formula  being 


udv  =  uv  -  vdu 


(2.40) 
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With  u  =  y  and  dv  =  (m  +y  )  ydy,  the  required  integration^  is 
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Evaluation  of  the  limits  gives 
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(2.37c) 


The  single  integral  of  this  result  can  be  further  integrated  to  give 
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At  this  point,  there  remain  the  integrals  that  arise  from  use  of 
the  parts  formula.  Fortunately,  however,  their  solutions  can  he  obtained 
from  the  integrals.  The  similarity  between  equation  (2.37d)  and  the 

previous  integrals  is  demonstrated  by  redesignating  the  summations 
variable  as  j  =  n+2.  The  relationship  becomes 


The  integrals  then  in  terms  of  J  are,  except  for  the  coefficients, 
identical  to  those  in  equations  (2.29a)  and  (2.35a)  for  ^c^4’ 

respectively.  This  condition  means  that  the  series  of  terms  to  be 
integrated  in  equation  (3*37)  is  the  same  as  the  series  beginning 

with  its  third  term.  Since  equation  (2.37e)  is  to  be  integrated  for 
n  =  3  to  n  =  n,  the  result  can  be  obtained  by  using  the  series 

from  its  5th  term  to  its  n+2nd  term.  This  series,  then,  with  the 
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coefficients  of  the  terms  modified  as  hy  equation  (3.36e)  gives  the 
required  integral.  The  results  are  given  in  the  following  summary  of 
the  y  and  z  integrations. 

Summary  of  First  Two  Integrations 

After  the  first  two  integrations,  the  relationship  describing 
the  function  I  can  be  expressed  as 

1=  (I,  +I„  +  I.  +  !,,)  +  (I-,  +  ^  _)  (2.4l) 
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Integration  With  Respect  to  Y 


At  this  point  in  the  development,  a  numerical  solution  must  he 
sought  since  the  integrand  of  equation  (2.19)  does  not  readily  lend 
itself  to  further  analytical  operations.  Also,  it  becomes  desirable 
to  solve  the  problem  on  a  digital  computer  and  to  make  nrunerical  evalua¬ 
tions  for  the  final  integration  variables. 

The  numerical  integrations  were  accomplished  by  applying  an 

appropriately  selected  member  of  the  closed  type  Newton-Cotes  quadrature 

33 

formulae.  These  formulae,  under  the  proper  conditions,  represent  quite 
accurately  the  integral  of  a  function  with  respect  to  one  of  its  varia¬ 
bles  provided  the  value  of  the  function  is  known  for  several  equally 


spaced  values  of  the  integration  variable.  The  specific  formula  employed 
for  the  numerical  computations  was 


■>  2 

J  ydx  =  (kly^  +  2l6y^  +  27y^  +  272y  +  27y; 


+  2l6yj.  +  4lyg)  -  2^ 


(2.42) 


which  was  found  to  give  accurate  results  with  an  acceptable  amount  of 
numerical  work.  In  this  equation,  y,  symbolizes  a  function  to  be  inte¬ 
grated  and  the  subscripted  y  variables  represent  the  ordinates  of  the 
function  at  seven  uniformly  spaced  values  of  the  independent  variable, 
indicated  by  x.  The  interval  between  ordinates  is  h.  The  symbols,  x 
and  h,  in  this  particular  problem  represent  Y  and  respectively. 

The  quantity,  9y^^^h^/l400,  is  the  error  incurred  by  use  of  the  formula, 

/  o  \ 

where  y^  is  the  eighth  derivative  of  the  function  with  respect  to  the 
variable  of  integration.  If  the  interval  is  small  (h^,  small)  and  the 
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function  is  not  exceedingly  irregular  (y'  known  everywhere  to  he 
small)  the  integral  for  most  practical  purposes  may  he  considered 
exact. These  conditions  were  met  except  for  r  and  Y  simultaneously 
approaching  R  and  zero,  respectively.  This  case  is  treated  separately. 

To  execute  the  numerical  calculations,  a  set  of  values  for  the 
radiation  system  must  he  chosen.  Since  the  theoretical  work  was  in¬ 
tended  to  he  complemented  hy  an  experimental  investigation,  initial 
calculations  were  made  for  variables  applying  to  the  experimental  sys¬ 
tem.  Specifically,  these  were  a  tube  diameter  of  0.043  ft.,  a  length 
of  1.04  ft.,  and  appropriately  selected  k  (attenuation  coefficient) 
values  of  0.5,  1.0,  4.0,  and  5.0  ft"^,  these  being  of  primary  interest 
to  the  experimental  study.  Later,  these  constants  were  varied  to  study 
their  effect  on  the  final  solution. 

With  the  above  constants  assigned,  the  development  may  he  con¬ 
tinued  hy  treating  the  terms  of  equation  (2.4l),  according  to  the 


(2.43c) 


where  the  integration  for  Y  is  defined  hy  equation  (2.19).  Inspection 
of  the  equation  shows  that  these  separate  integrals  have  special  signi¬ 
ficance.  is  the  collection  of  terms  not  having  a  k-coefficient  andj 
hence,  represents  a  pseudo  view-factor  between  the  radiating  wall  and 
the  absorbing  particles.  Physically,  the  result  from  the  treatment  of 
may  be  envisioned  as  an  evaluation  of  the  total  direct  transfer  of 
heat  with  no  particle  shading.  This  function  also  has  the  characteris¬ 
tic  of  being  undefined  for  m  =  0;  it  requires  special  treatment..  The 
second  integral,  I  ,  represents  the  effect  of  radiant  energy  attenuation 

Zj 

by  particle  shading.  Since  attenuation  reduces  the  exposure  of  the 
particle  cloud,  is  always  negative  and  subtracts  from  J^.  Finally, 

I  is  defined  for  all  values  of  m  and  its  numerical  evaluation. is  much 

Zj 

different  for  computer  programming  than  it  was  for 

Numerical  calculations  were  initiated  with  the  evaluation  of  the 
attenuation  function,  J^,  which,  as  previously  denoted,  was  accomplished 
by  use  of  equation  (2.42).  The  procedure  using  the  integral  formula  to 
establish  the  final  result  was  to  calculate  ordinates  (l^)  at  seven 
uniformly  spaced  values  of  Y  within  the  range  of  Y  across  which  the  in¬ 
tegration  was  being  performed  and  to  incorporate  these  ordinates  in  the 
integral  formula.  The  range  of  values  of  Y  across  which  the  formula 
was  applied  was  variable,  and  the  procedure  could  be  repeated  as  many 
times  as  required  to  yield  the  accuracy  desired,  i.e. ,  the  increment  h 
(or  ay)  of  the  error  function  was  arbitrary.  For  this  investigation,  • 
the  integral  formula  was  applied  twice,  once  in  each  interval  ’ 

0  ^  Y  ^  tt/2  and  n/2  ^  Y  n,  and  the  total  integral  was  the  sum  of  the 


results  for  each  interval. 


Calculated  results  for  applying  to  the  experimental  system  are 
presented  in  Figures  7  through  10.  The  data  are  given  for  radii  of  0.0, 
0.25R,  0.50R,  0.75R,  and  l.OOR  since  results  for  these  points  were 
required  for  the  final  integration  with  respect  to  radius.  The  func¬ 
tional  behavior  of  the  curves  shown  on  the  figures  was  employed  in 
studying  the  correctness  of  the  first  two  integrations  (y  and  z)  which 
resulted  in  a  lengthy,  complex  function.  Physically,  the  curves 
represent  the  attenuation  within  an  elemental  volume  receiving  radiation 
from  a  wall  element  where  the  linear  dimension  of  each  element  along  the 
axis  of  the  cylinder  has  been  integrated;  the  magnitude  of  the  ordinates 
are  measures  of  the  associated  reductions  in  radiation  intensity  of  rays 
passing  through  the  absorbing  element.  The  numerical  values,  although 
being  for  the  experimental  system,  are  generalized  results,  since  they  are 
functions  of  L/R,  the  attenuation  coefficient  k,  and  the  cylinder  radius 
R.  The  magnitude  of  the  parameter  kR  (attenuation  coefficient  times  the 
radius)  relative  to  unity  affects  the  appearance  of  the  curves,  however. 
The  influence  of  the  parameter  is  subsequently  discussed.  The  results 
presently  considered  are  for  small  kR  values  of  0.1  and  less.  As  shown 
by  the  figures,  there  is  no  angular  dependence  for  r  =  0,  as  would  be 
expected.  At  r  =  R,  on  the  other  hand,  the  attenuation  is  greatest  when 
the  two  differential  elements  are  diametrically  opposed  (f  =  n)  and  mono- 
tonically  decreases  with  f  until  the  element  occupies  the  same  position 
(f  =  O).  At  that  point,  there  is  no  attenuation  within  the  differential 
volume.  At  intermediate  radii,  the  attenuation,  conversely,  is  largest 
at  Y  =  0,  and  is  a  minimum  at  some  angle,  The  angle  at  which  the 

minima  occiir  is  a  function  of  radius  and  attenuation  coefficient.  The 
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Figure  7*  Ij’  as  Function  of  ip  and  Radial  Position  for  k=  0.5* 
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Figure  8.  ly  as 
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Figure  9*  Ij;  Fixnction  of  ip  and  Radial  Position  for  k=  l+.O 


Figure  10.  as  Function  of  ip  and  Radial  Position  for  k=  5*0 


loci  of  these  points  are  indicated  on  the  figures  by  dotted  lines. 

That  the  attenuation  should  be  largest  at  the  point  of  nearest 

approach  might,  at  first,  seem  illogical.  This  is  correct,  however, 

for  attenuation  is  a  first  order  process,  as  evidenced  by  the  applica- 

kx 

bility  of  Beer's  law,  i.e.,  i^^  =  i^e"  .  Hence  the  decrement  of  radia¬ 
tion  within  a  differential  volume  is  proportional  to  the  intensity 
(di  =  -ki  dx).  Since  the  intensity  of  a  beam  is  largest  at  Y  =  0, 

'  X  X 

the  attenuation  there  is  a  maximum  for  all  radii.  This  maximum  for  the 
r  =  R  curve  is  -rrkL,  but,  because  of  reasons  given  below,  this  particular 
curve  must  pass  through  zero  first  and  the  -rrkL  limit  has  no  significance. 

The  skewed,  sinusoidal  behavior  of  the  1^,  curves  for  constant 
radius  is  due  to  Lambert's  cosine  law.  Hence,  when  cos3^  decreases 
the  radiant  intensity  of  a  beam  must  be  reduced  and,  because  of  first- 
order  behavior,  the  attenuation  is  decreased.  The  minimum  attenuation 
then  occurs  where  cos0^  is  the  smallest.  This  position  is  located  where 
the  radiant  beam  is  tangent  to  the  radial  circle  being  considered.  On 
either  side  of  this  value  the  angle  3^  must  decrease  thus  giving  a 
larger  intensity  (i^^^  =  i^^^^  cos3^).  The  loci  of  these  minima  are 
given  by 

r  =  cos-^(|) 

m  ^ 

The  angular  function  described  by  equation  (2.44),  however,  is 
not  the  only  factor  in  determining  the  position  of  minimum  attenuation. 

It  may  be  seen  from  the  figures  that  the  actual  minimum  attenuation 
occurs  at  a  larger  angle  than  the  predicted  by  equation  (2.44).  For 
example,  when  k  is  -5-0,  and  r  =  0.50R,  the  minimum  should  be  at  tt/3, 


but  is  actually  at  5tt/12.  Further  inspection  of  the  curves  show  that 
the  discrepancy  is  a  function  of  k  and  r.  At  large  r  and  small  k 
values,  equation  (2.44)  approximates  the  correct  result  very  closely; 
as  r  approaches  zero,  however,  the  equation  does  not  hold  at  all  as 
the  minimum  attenuation  point  approaches  n.  The  shift  in  the  angle 
giving  smallest  attenuation  to  values  larger  than  predicted  by 
equation  (2,44)  is,  due  to  the  first  order  dependence  of  attenuation 
on  the  distance  between  the  emitting  wall  element  and  the  differential 
volume  of  absorbing  aerosol.  Equation  (2.44),  then,  is  exact  for  con¬ 
stant  intensity  but  it  must  be  amended  for  decreasing  intensity  with 
distance.  The  exact  effects  of  changing  intensity  are  not  readily 
calculable;  it  is  apparent,  however,  that  the  angular  dependence  is 
dominant  for  large  radii  (r  R)  while  for  small  radii  (r  O)  the 
distance  of  separation  is  more  important.  Each  of  the  factors  are 
completely  controlling  for  r  =  R  and  r  =  0,  respectively.  Special 
attention  should  be  given  to  the  fact  that  at  r  =  R  the  angle  with 
minimum  attenuation  is  exactly  zero  and  that  the  curve  steadily 
decreases  with  Y  and  approaches  zero  as  a  limit*  At  Y  =  0  and  r  =  R, 
there  is  also  a  discontinuous  limit  of  -nkL,  but  this  point  is  of  no 
consequence  in  evaluating  the  area  between  the  r  =  R  curve  and  the 
axis.  The  -rrkL  point  is  the  limit  that  all  other  radial  curves 
approach  with  increasing  r  and  decreasing  Y. 

The  distance  between  the  two  elemental  areas  exchanging  radiation 
influences  the  above  particular  results  for  Figures  7  through  10  slight¬ 
ly,  since  the  parameter  kR  is  small.  The  curves  are  characteristically 
determined  by  the  angular  function  given  by  equation  (2.44),  When  kR 
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is  large,  however,  attenuation,  as  it  is  affected  by  the  distance  of 
separation,  dominates  the  result  and  its  dependence  on  is  almost 
completely  obscured.  Some  residual  dependence  on  the  angle  remains 
near  the  radiating  wall  element.  Supplementary  calculations  are  pre¬ 
sented  in  Figures  11  and  12  to  emphasize  these  characteristics.  There 
it  may  be  observed  that  for  small  values  of  Y  the  curves  behave  simi- 
lax  to  those  already  described  for  small  kR.  As  Y  or  r  changes  so  as 
to  separate  further  the  radiating  wall  element  and  the  absorbing 
volume  element,  the  energy  attenuated  by  the  latter  is  drastically 
reduced  because  of  the  rapidly  decreasing  intensity  therein  received. 
Hence,  the  curves  for  any  radius  approach  zero  with  increasing  t  . 

It  will  also  be  noted  that  the  position  of  the  r  =  0  curve  is  rela¬ 
tively  much  closer  to  the  axis  for  very  large  kR  values.  The 
minimum  ordinate  for  all  curves  (except  r-^)  is  at  Y  =  tt. 

When  the  I  functions  are  integrated  with  respect  to  Y  the 
result  again  has  physical  significance  and  may  be  thought  of  as  being 
the  energy  attenuated  from  radiation  passing  through  a  cylindrical 
shell  with  length  L,  radius  r,  and  thickness  dr.  The  calculated 
results  designated  as  are  given  in  Figure  13.,  There  it  may  be 
observed  from  the  magnitude  of  the  ordinates  that  shells  near  the  center 
of  the  tube  screen  more  energy  per  unit  circiamferential  length  than 
shells  near  the  wall.  This  result  must  be  interpreted  on  a  per  unit 
length  basis  since  the  variation  of  r  (or  circmference)  is  yet  to  be 
considered  (next  integration).  Again  the  numerical  data  for  the  J^,  cal¬ 
culations  really  pertain  to  the  experimental  system  but  were  employed  as 
generalized  functions  of  L/R  and  k  values,  and  radial  position. 
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Figure  11.  as  Function  of  ijj  and  Radial  Position  for  kR=  2.0. 
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Figure  12.  as  Function  of  ip  and  Radial  Position  for  kR=  1.0 


Figure  13.  J  as  Function  of  ip  and  Radial  Position  and  Attenuation 
Coefficients. 


The  evaluation  of  the  integral  is  more  involved  than  be¬ 
cause  the  function  becomes  undefined  for  approaching  zero  and  the 
numerical  procedure  using  equation  (2.42)  becomes  invalid.  The  situa¬ 
tion  can  be  overcome,  to  some  extent,  by  reducing  the  AY  interval  and 
using  more  applications  of  the  integration  formula.  Ultimately,  how¬ 
ever,  the  limiting  case,  where  m  vanishes,  requires  a  complementary 
analytical  analysis.  The  characteristics  of  the  integrand  (l^) 
curves  are  demonstrated  in  Figure  l4.  They  are  functions  of  L/R  and 
R  only.  On  the  figure  the  curves  may  be  noted  as  being  elliptic  with 
the  centerline;  the  r  =  R  curve  is  a  special,  limiting  case.  While 
the  data  shown  are  for  a  particular  (L/r)  of  50,  it  may  be  demonstrated 
from  equation  (2,42)  that  the  curves  are  functions  of  Y  and  (l/r)  and 
the  configurations  given  are  general.  For  smaller  (l/r)  ratios,  however, 
the  ordinates  for  small  Y  are  not  quite  as  large  and  the  nearly  horizon¬ 
tal  portions  of  the  curves  (r  ^  O)  are  less  pronounced.  The  function 
for  r  =  R  is  undefined  at  Y=0;  and  at  Y  =  rr,  it  has  the  value 

C/(i)%  4 2] . 

The  procedure  for  evaluating  the  integral  of  with  respect  to 
Y  was  to  begin  operations  with  r  0  and  to  modify  the  calculations  for 
handling  the  erratic  behavior  with  increasing  radius.  At  r  =  0  the 
integral,  equation  (2.43a),  readily  reduces  to 

•^l,r.O  =  E  I  if 

o 

and 

(2.k6) 
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Figure  lU.  8.S  Function  of*  ^  nnd  RsKiia.!  Posit-ion* 


All  other  radial  positions  require  numerical  treatment  and,  as 
before,  for  the  calculation  of  J^,  the  Newton-Cotes  numerical  integra¬ 
tion  formula,  equation  (2.42),  was  employed.  For  the  range  of  radii 
0  ^  r  ^  0.6R,  the  integral  formula  was  applied  eight  times  to  cover 
the  angular  range  of  Y  from  zero  to  tt,  and  calculations  of  were 
made  for  radii  of  O.IR,  0.2R,  0.3R,  0.4R,  0.5R,  and  0.6R.  At  still 
larger  radii,  calculations  of  were  made  for  equally  spaced  radial 
values  of  0.66R,  0.73R,  0.86R,  and  0.93R.  The  integral  formula  was 
used  here  eight  times  in  each  of  the  angular  intervals  0  ^  Y  ^  0.05, 
0.05  ^  ’i'  ^  0.1,  and  0.1  s  Y  ^  n. 

Finally,  because  the  integrand  increases  without  limit  as  r 
approaches  R,  the  numerical  calculations  have  to  be  supplemented  with 
an  analytical  development.  The  argument  is  begun  by  considering  the 
integral  in  equation  (2.43)  to  be  con^josed  of  three  parts. 
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with  r  approaching,  but  not  equal  to,  R  and  taken  to  be  small.  For 
these  conditions,  the  result  can  quite  accurately  be  approximated  by 


57 


In  the  limit  this  expression  is  exact.  Also  for  small  angles,  the 
cosine  can  be  taken  as 


cos  Y  =  1  -  0.5Y 


2 


(2.49) 


•which  is  the  first  two  terms  of  the  power  series  of  cos  Y.  It,  too, 
becomes  exact  in  the  limit.  With  these  approximations  the  integral  when 
r  is  approaching  R  can  be  e-valuated: 
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or 


Jl  =  (2)(|)(^^)tan"^  (|^)  -  (f)(^=)  Sinh’^  (|^) 
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The  following  limits  are  then  considered: 
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(2.55) 


The  value  of  for  r=R  is,  finally 


n 


2  .2x2. 


dY 


(2.56) 


Since  the  integral  of  the  later  equation  does  not  become  undefined  for 
any  value  of  r  or  Y,  it  can  be  evaluated  numerically.  The  r=R  asymptote 
of  Figure  l4  refers  to  the  experimental  system  to  be  described  subse¬ 


quently. 

The  complete  results  for  and  given  in  Figure  15  for  all 
values  of  r.  The  function  represented  in  this  case  is  the  total  heat 
absorption  per  unit  circvimferential  length  occurring  within  a  cylindri¬ 
cal  shell  of  radius  r,  length  L,  and  thickness  dr.  It  may  be  noted  that 
the  absorption  of  energy  for  a  large  L/R-ratio  tube  is  almost  indepen¬ 
dent  of  radius  whereas  the  attenuation  as  shown  in  Figure  13  has  very 
pronounced  variation.  The  characteristics  of  each  function,  and  J^, 
are  probably  strong  functions  of  the  parameter  kR  and  L/R  ratio,  however. 
Integration  with  respect  to  r. 

The  final  integral  with  respect  to  r,  according  to  equations 
(2.19)  and  (2.43)  is 
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Expanded  J,  Sea 


(2.57) 


R 

=  2TTi,  CS  R  f(j,+j„)rdr 
bnw  V  J  ^  1  Z 
o 

and 


=  2nL  CS  R(F'+f;) 

bnw  V  1  Z 


(2.58) 


where  is  the  total  direct  heat  transfer  to  a  cloud  of  particles. 

Like  the  immediately  preceding  integration  with  respect  to  Y,  this  one 
was  also  performed  numerically  and  the  and  functions  were  evalu¬ 
ated  at  such  radial  points  that  a  closed  Newton-Cotes  quadrature  for¬ 
mula  again  could  he  employed.  The  function  was  integrated  by  apply¬ 
ing  the  seven  point  formula  [equation  (2.42)],  once  in  the  interval 
0  r  s  0.6R  and  once  in  the  interval  0.6R  ^  r  s  l.OOR.  Because  of  the 
much  longer  computer  time  involved  in  calculations  leading  to  less 
radial  points  were  computed  and  a  five  point  integration  formula  was 
employed  for  the  evaluation  of  F^.  The  five-point,  closed- form, 
Newton-Cotes  equation  is 

^2  (6)  7 
J  ^^^1  ^  ^^^2  %4'^~~  (2.59) 

where  h  is  the  increment  for  the  integration  x,  y  is  the  function  to 
be  integrated,  the  subscripted  y's  are  the  ordinate  values  required, 
and  [8y(^^h'^]/945  is  the  error  with  y^^^  being  the  sixth  derivative  of 
the  integrand. 

The  final  results  for  are  presented  in  the  discussion  of 
results  (Chapter  IV)  where  the  calculated  heat  transfer  is  compared  with 
experimental  data  and  generalized  calculations  are  evaluated. 
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CHAPTER  III 


EXPERIMENTAL  INVESTIGATIONS 

Description  of  Apparatus 
Procedure  and  General  Operation 

Previous  studies  of  heat  transfer  to  particle  clouds  flowing  in 

cylindrical  enclosures  have^  in  general j  emphasized  the  conduction  and 

convection  processes ,  and  the  contributions  of  radiation ^  when  signifi- 

35  3^  3V 

cant,  have  been  estimated  by  theoretical  calculations  alone.  ’  ’ 

Experimental  studies  considering  dispersed  systems  usually  treat  luminous 
flames  for  which  the  desired  radiation  properties  are  derived  from  opti¬ 
cal  pyrometer  measurements.  Hence  the  direct  absorption  of  radiation  by 
particle-gas  systems  has  not  been  determined.  In  this  investigation, 
therefore,  measurements  were  made  to  provide  precise  data  on  the  trans¬ 
fer  of  radiant  energy  to  systems  consisting  of  nonluminous  particles 
well- characterized  as  to  particle  size,  surface  area,  and  concentration 
suspended  in  a  gas.  The  first  problem  was  to  develop  an  experimental 
system  in  which  a  particle  cloud  could  be  exposed  to  a  radiant  field 
having  conduction  and  convective  heat  transfer  reduced  to  a  minimum  and 
the  energy  could  be  absorbed  by  the  particle  cloud  in  such-  a  way  that 
its  magnitude  might  be  determined  by  enthalpy  balances  before  and  after 
exposure. 

The  final  experimental  system  exposed  a  particle  cloud  to  a 
radiant  heat  flow  by  passing  it  through  a  quartz  tube  which  transmitted 
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to  its  interior  a  significant  portion  of  the  radiation  falling  on  it 
externally.  The  quartz  tube  was  protected  from  extraneous  heat 
transfer  effects  other  than  radiation  by  enclosing  it  in  an  evacuated 
space,  the  outer  surface  of  which  was  air-cooled.  The  surrounding, 
protective  tube  was  also  quartz.  The  particle  cloud  was  generated  at 
a  constant  rate  and  pre-conditioned  to  a  constant  known  temperature. 
Then,  after  passing  through  the  radiation  field,  the  cloud  was  cooled 
calorimetrically  to  determine  the  heat  absorbed.  A  filter  mechanism 
removed  the  suspended  particles  from  the  air  and  returned  them  to  an 
aerosol  generator  for  re-use.  The  apparatus  was  thus  designed  for 
continuous  operation.  Repeated  experiments  could  be  performed  with  a 
limited  inventory  of  powder  (one-half  pint)  which,  in  some  cases, 
required  lengthy  preparations  meaning  that  it  could  not  readily  be 
obtained  in  large  quantities.  Figure  l6  is  a ‘phbtograph  of  the 
assembly,  and  Figure  17  is  a  diagramatic  representation  showing  the 
functions  of  the  principal  components. 

Description  of  Principal  Components 

Particle  Cloud  Preparation.  The  particle  cloud  or  aerosol  was 
generated  with  the  unit  identified  as  an  aerosol  generator  in  Figure 
17  and  drawn  in  detail  in  Figure  I8.  As  shown  in  section  A-A  of  the 
latter  drawing,  a  powder  charged  to  the  generator  collected  behind  a 
retaining  partition  in  a  slowly  rotating  container.  A  thin  layer  of 
powder  was  thus  produced  along  the  bottom  of  the  container  due  to  the 
latter *s  moving  continuously  beneath  the  retaining  partition.  This 
layer  was  carried  by  the  revolving  surface  under  nozzle  pickups  where 
the  powder  was  entrained  in  air  passing  upward.  Particle  deagglomera- 


Figvire  IT.  Diagramatic  Representation  of  Heat  Transfer  Apparatus. 
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Figure  l8.  Particle  Cloud  Generator, 


tion  resulted  from  the  shear  forces  generated  hy  the  flow  through  the 
small  passages;  the  pressure  drop  across  them  was  approximately  10 
psig.  The  newly  formed  particle  cloud  next  passed  through  a  length  of 
pipe  (L/D=100)  to  allow  the  streams  from  the  nozzle  to  mix  thoroughly 
and  to  reach  a  steady  flow  condition.  The  entrance  length  was  enclosed 
in  an  air  cooled,  constant  temperature  jacket  to  bring  the  aerosol  to 
a  known  temperature  for  use  later  in  computing  an  overall  heat  balance 
on  the  suspension  system. 

Numerous  exploratory  experiments  were  required  to  achieve  a 
suitable  de agglomerating  nozzle  to  handle  a  large  flow  rate  and  still 
produce  a  uniformly  dispersed  cloud  with  steady  flow  characteristics. 
Originally,  a  venturi  reducing  the  channel  diameter  to  about  0.03  inch 
was  employed  at  the  point  the  powder-air  mixture  issued  from  the 
generator  assembly.  At  that  time  the  powder  was  made  airborne  by  en¬ 
trainment  with  air  passing  through  a  perforated  plate  that  covered  the 
inlet  to  the  aerosol  conduit.  In  practically  every  case  heavy  agglo¬ 
merates  would  form  a  fluidized  bed  in  the  pipe  leading  to  the  venturi. 

The  venturi  was  then  placed  inside  the  aerosol  generator  to  avoid  the 
fluidizing  effects,  but  now  the  powder  entrainment  was  unsteady.  The 
most  satisfactory  arrangement  employed  multiple,  small-diameter  passages 
inside  the  generator  which  served  both  as  pickup  tubes  and  as  deagglo¬ 
meration  nozzles.  Suspensions  up  to  0.05  cent  solids  by  volume  -- 
equivalent  to  75  per  cent  by  weight  for  dense  materials  --  were  produced. 

The  aerosol  conduit  was  0.5  inch,  in  diameter,  A  larger  system 


could  not  be  operated  without  overloading  the  compressed  air  facility 
sixpp lying  the  laboratory. 
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Radiant  Heating  Section.  After  preparation,  the  particle  cloud 
was  e^cposed  to  a  radiant  field  generated  by  Kanthal  wire,  (Kanthal 
Company,  Stanford,  Connecticut)  electric  resistance  heaters  supported 
in  grooves  spiraling  up  the  inner  walls  of  a  refractory  enclosure  con¬ 
structed  with  Johns -Mans ville  23  refractory  brick.  Two  heaters  were 
employed,  one  for  the  upper  and  the  other  for  the  lower  half  of  the 

furnace.  Each  was  composed  of  30  feet  of  0.333  ohm  per  foot  wire 

wound  into  a  0.375  inch  diameter  coil.  The  coil  helices  were  spaced 
closer  to  one  another  at  the  ends  of  the  furnace  to  compensate  for  the 
heat  losses  there  and  thus  to  give  a  more  uniform  temperature  over  the 
heated  area.  Further,  the  temperatures  of  the  two  coils  were  controlled 
independently;  this  arrangement  reduced  the  temperature  differences 
resulting  from  convective  currents  within  the  furnace.  A  maximum  of 

k  kw  could  be  supplied  to  the  furnace  and  2300^  F  was  the  maximum 

operating  temperature  for  the  resistance  wire.  The  directly  irradiated 
aerosol  path  was  12.5  inches  in  length.  Figure  19  per sent s  a  drawing 
of  the  assembly. 

The  aerosol  conduit  passing  through  the  radiant  field  was  a 
13  mm  ID  quartz  tube  having  a  1.0  mm  wall  thickness.  Surrounding  this 
tube  were  two  others,  the  inner  one  25  imn  ID  with  a  1.0  mm  wall  thick¬ 
ness  and  the  outer  one  mm  ID  with  a  1.5  nun  wall.  The  middle  tube 

was  sealed  to  the  aerosol  conduit  and  the  annular  enclosure  thus  formed 
was  evacuated  to  a  pressure  of  approximately  1.0  micron  of  mercury  to 
reduce  heat  transfer  by  conduction  and  convection  across  the  enclosed 
space  into  the  suspension.  Air  was  circulated  in  the  outer  annulus  to 

68 


protect  the  outer  quartz  tubes  from  the  high  furnace  temperatures.  The 
air  temperature  in  the  outer  tube  could  be  held  below  200  F  with  air 
supplied  from  a  suction  fan.  Because  of  the  vacuum  between  this  enclo¬ 
sure  and  the  aerosol  conduit,  very  little  heat  could  be  transferred 
other  than  by  radiation  to  the  aerosol  system  being  studied.  Some  heat 
could  be  transferred,  of  course,  where  the  vacuum  jacket  joined  the 
aerosol  conduit,  and  some  energy  was  absorbed  by  the  quartz  tube  carry¬ 
ing  the  aerosol.  Nevertheless,  these  heat  losses  were  small,  and 
experimental  corrections  could  be  applied  to  evaluate  the  radiation 
transferred  to  the  particles. 

Since  the  radiation  had  to  pass  through  three  quartz  walls  having 
a  total  thickness  of  3*5  millimeters  before  reaching  the  aerosol,  the 
radiation  was  considerably  altered.  An  estimation  of  this  change  was 
made  from  other  experiments  using  the  same  type  quartz  as  employed  in 
the  furnace  section.  Figure  20  presents  data  for  a  2  millimeter  thick- 

oQ 

ness  as  supplied  by  the  manufacturer.  Reflection  losses  were  not 
included,  but  from  other  measurements  they  were  found  to  be  small. 

As  shown,  the  transmission  efficiency  was  well  over  90  per  cent  in  the 
0.3  to  3.5  micron  wave  length  interval  except  for  a  narrow  dip  at  2.7 
microns,  and  the  transmission  was  nil  over  4.7  microns.  Most  of  the 
radiation  exposure  for  the  aerosol,  then,  was  in  the  range  between 
0.3  and  3.7.  Since  the  temperature  of  the  furnace  was  near  2500^  R, 
the  shorter  wave  lengths  were  further  limited  to  those  above  about 
1  micron.  The  peak  intensity  was  around  2  microns.  The  effective 
radiation  field,  then,  was  almost  completely  between  1.0  and  3.0 
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microns.  The  average  intensity  of  this  field  was  determined  by 
ejqjeriment . 

Heat  Exchange  Section.  The  purpose  of  the  experimental  work  was 
to  measure  the  heat  gained  by  the  aerosol  system  from  the  radiation 
within  the  furnace.  It  would  have  been  desirable  to  make  the  measure¬ 
ments  directly  as  the  aerosol  passed  from  the  furnace,  but  experiments 
of  this  type  did  not  give  satisfactory  results.  Heat  was  conducted  as 
well  as  radiated  into  the  downstream  conduit,  and  this  was  found  to  lead 
to  ambient  effects  that  obscured  the  desired  measurements.  Equilibrium 
was  both  slowly  attained  and  unstable.  Changes  in  temperature  could  be 
detected  when  particles  were  added  to  the  system,  but  interruption  of 
the  powder  flow  frequently  did  not  return  the  system  to  its  original 
thermal  state. 

Better  results  were  obtained  when  the  aerosol  conduit  was  incor¬ 
porated  in  a  concentric,  single  pass,  tube -and- shell  heat  exchanger. 

The  latter  unit,  which  was  well  insulated  from  the  surroundings,  served 
as  a  controlled  heat  sink  between  the  furnace  and  a  downstream  point 
where  the  aerosol  temperature  was  measured.  The  steady  removal  of  heat 
stabilized  the  thermal  behavior  of  the  apparatus  because  disturbances 
were  readily  absorbed  by  the  heat  exchanger  coolant  and  were  not  per¬ 
mitted  to  remain  as  sensible  heat  in  the  apparatus  to  be  removed  by  the 
aerosol  stream.  Hence,  measurements  at  a  given  time  were  more  precise 
and  were  not  dominated  by  transient  thermal  conditions.  A  sectional 
drawing  of  the  heat  exchanger  is  shown  in  Figure  21.  The  enclosed  con¬ 
duit  conveying  the  particle-gas  suspension  was  0.5  inch  in  diameter  and 
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four  feet  long.  The  conduit  was  surrounded  hy  a  one-inch  tube  through 
which  a  cooling  medium  flowed  countercurrent  to  the  suspension.  Either 
or  water  could  be  the  coolant ^  and  a  choice  could  be  made  on  the 
basis  of  which  had  the  appropriate  heat  capacity  to  give  a  readily  mea¬ 
sured  temperature  rise  in  passing  through  the  system.  For  the  parti¬ 
cular  suspensions  studied,  air  gave  better  results.  A  vacuum  chamber 
surrounded  the  heat  exchanger  to  insulate  it  from  the  ambient  room  con¬ 
ditions.  The  Inside  wall  of  this  chamber  was  lined  with  aluminum  foil 
and  the  outer  surface  of  the  heat  exchanger  shell  was  polished  to 
reduce  radiation  losses.  End  pieces  were  made  as  thin  as  possible  to 
minimize  the  losses  there,  and  insulation  was  provided  to  stabilize 
the  temperatures  to  which  the  highly  conducting  end  pieces  were 
exposed.  The  assumption  was  made  that  the  localized  heat  disturbances 
at  the  ends  were  approximately  constant  during  an  experiment  and  could 
be  eliminated  by  heat  balance  calculations. 

Measurements  of  the  inlet  and  exit  temperatures  of  the  cooling 
medium  and  the  aerosol  exit  temperature  were  made  with  copper-constan- 
tan  thermocouples  arranged  as  shown  on  the  figure.  These  temperatures 
together  with  the  mass  flow  rates  of  the  two  streams  permitted  making 
a  heat  balance  for  each  and,  from  these,  evaluating  the  heat  content  of 
the  aerosol.  The  temperature  of  the  aerosol  itself  was  also  measured 
by  a  thermocouple  installed  6.5.  inches  from  the  entrance  to  the  heat 
exchanger.  By  assuming  a  linear  temperature  gradient  for  the  aerosol 
along  the  heat  exchanger,  this  measurement  could  be  corrected  to  give 
the  aerosol  temperature  at  the  entrance.  A  supplementary  heat  balance 
made  exclusively  on  the  aerosol  system  was  then  possible.  A  radiation 
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shield  in  the  aerosol  streajn  prior  to  its  entrance  to  the  heat  exchanger 
stopped  downstream  radiation  in  the  heat  exchanger  and, because  of  the 
turbulence  generated,  induced  a  more  rapid  adjustment  of  particle  and 
gas  temperatures.  This  enhanced  the  reliability  of  the  lower  tempera¬ 
ture  measurement.  The  shield  consisted  of  a  round-edge  flow  nozzle 
gently  reducing  the  conduit  cross-section  to  about  one  fourth  its  ori¬ 
ginal  size  and  a  downstream  conical  section  with  its  apex  at  the  dis¬ 
charge  side  of  the  orifice  to  cut  off  direct  longitudinal  transmission 
of  radiation. 

Particle  Concentration  Monitoring  and  Sampling  Methods.  Upon 
leaving  the  heat  exchanger,  the  particle  cloud  passed  through  a  return 
bend  at  the  top  of  the  apparatus.  In  the  downward  path  it  was  moni¬ 
tored  for  consistency  of  concentration  by  means  of  light  transmission 
measurements  made  as  shown  in  Figure  17  (page  65  )  The  principal  com¬ 
ponent  of  the  device  was  the  special  electronic  circuit  diagrammed  in 
Figure  22;  it  provided  a  constant  voltage  for  the  lamp  which  generated 
a  light  beam  of  substantially  constant  intensity.  The  lajnp  was  posi¬ 
tioned  so  that  it  illuminated  a  photoelectric  cell  diametrically  across 
the  aerosol  path.  Since  the  light  transmitted  across  this  space  was  a 
function  of  the  particle  concentration,  the  photocell  current  was  a 
measure  of  the  number  of  particles  present. 

Light  transmission  measurements  were  intended  to  serve  primarily 
as  a  means  of  monitoring  the  fluctuations  in  particle  concentration 
during  the  course  of  experiments.  In  theory,  the  electronic  instrument 
could  have  been  calibrated  for  each  powder  and  the  transmissivities 
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could  have  been  then  converted  directly  into  concentration.  In  practice 
however,  considerable  uncertainty  is  introduced  by  this  approach.  For 
instance,  the  illujnination  varies  with  the  cleanness  of  the  conduit 
windows.  The  device,  therefore,  was  not  reliable  for  absolute  measure¬ 
ments,  but  it  was  useful  for  detecting  concentration  changes  during  an 
experiment  and  for  assuring  steady  operation  of  the  aerosol  generator. 

An  absolute  measure  of  the  particle  cloud  concentration  was 
obtained  by  direct  sampling.  The  total  powder  output  of  the  aerosol 
generator  was  measured  for  the  four  different  speeds  available  in  its 
operation.  These  results,  together  with  measurements  of  the  air  flow 
through  the  system,  then  permitted  the  particle  concentration  to  be 
calculated.  The  experimental  procedure  was  to  close  off  the  aerosol 
from  its  normal  path  around  the  heat  transfer  circuit  and  to  open 
simultaneously  a  sample  tube  leading  to  an  extraction  thimble  which 
collected  the  entrained  particles.  After  a  short  interval  of  time, 
usually  about  30  seconds,  the  aerosol  flow  was  diverted  back  to  its 
normal  path  and  the  sample  was  removed  and  weighed.  The  results  were 
corrected  to  account  for  the  powder  retained  in  the  sampling  connector. 
Between  6  and  10  determinations  were  made  for  each  concentration  level, 
and,  to  be  considered  acceptable  tests,  agreement  between  results  had 
to  be  within  ±  10  per  cent;  they  frequently  were  within  ±  5  per  cent. 

A  Staplex  industrial  air  sampler  (The  Staplex  Company,  777  Fifth 
Avenue,  Brooklyn  32,  New  York)  maintained  constant  pressure  across  the 
thimble  filter  to  minimize  the  interference  of  the  sampling  procedure 
with  the  normal  operation  of  the  system. 
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Powder  Collection  and  Return  Mechanism «  Finally  the  aerosol 
entered  the  sleeve  filter  (shown  in  Figure  1?)  where  the  powder  was 
separated  from  the  air  stream.  The  filter  medium  was  4  inches  in 
diameter  and  had  22  inches  of  useful  length.  The  filter  was  traversed 
outside  hy  a  circumferential  air  jet  which  dislodged  the  filtered 
material  and  provided  continuous  operation  of  the  filter.  The  traver¬ 
sing  action  was  obtained  with  a  guide  that  followed  a  right-hand  helical 
groove  along  a  cylinder  in  one  direction  and,  at  the  end,  reversed  to 
follow  a  left-hand  helix  to  the  opposite  end.  The  shifting  of  the  guide 
from  one  helix  to  the  other  was  accomplished  hy  a  connecting  groove  at 
each  end  of  the  cylinder.  The  guide  had  sufficient  length  to  avoid 
changing  helixes  where  they  cross.  The  mechanism  was  driven  hy  a  motor 
and  gear  speed-reducing  unit. 

Two  discharge  valves  helow  the  filter  transferred  the  collected 
powder  hack  to  the  aerosol  generator  for  reuse.  These  valves,  shown 
in  Figure  23^  were  arranged  in  series  and  were  operated  so  that  one 
valve  was  always  closed  to  prevent  escape  of  particles  from  the  system. 
The  valves  were  operated  manually. 

Controls  and  Instrumentation.  Several  control  and  measurement 
devices  were  required  for  the  operation  of  the  apparatus.  Figure  24 
is  a  photograph  of  them.  The  right  side  of  the  panel  contained  the 
necessary  rotameters  for  measuring  the  various  air  and  water  flow  rates. 
The  rotameter  at  the  extreme  right  measured  the  air  flow  rate  supplied 
to  the  heat  exchanger  for  cooling  the  particle  cloud,  while  the  inter¬ 
mediate  rotameters  metered  air  to  the  particle  cloud  generator  for 
producing  the  aerosol.  The  remaining  rotameter  was  to  meter  water  to 
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Figure  23.  Powder  Discharge  Mechanism. 
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Figure  24.  Instrument  and  Control  Panel. 


the  heat  exchanger  when  water  served  as  the  cooling  medinm.  The  valves 
beneath  each  rotajneter  controlled  the  flow  through  the  meter  above. 

The  left  part  of  the  panel  contained  the  power  and  electronic 
equipment.  Across  its  lower  portion  was  a  master  switch  for  the 
apparatus  and  the  two  variable  auto-transformers  that  supplied  power 
to  the  furnace  coils.  Above  was  the  switch  panel  for  the  various 
power  sources  throughout  the  apparatus.  And  at  the  right- center  was 
the  thermocouple  selection  switch  for  the  temperature  indicating  system. 
The  recorder  (Speedomax  W,  Leeds  and  Northrup,  Philadelphia,  Pa.)  indi¬ 
cated  the  furnace  temperature. 

Unmounted  equipment  included  a  model  MR  multi-range  and  a  model 
SR  Sargent  recorder  (E,  H.  Sargent  &  Co.,  Chicago  30,  Ill.),  a  Leeds 
and  Uorthrup  No.  8686  potentiometer,  and  a  Leeds  and  Northrup  No.  2430-A 
galvamometer.  These  units  measured  the  required  stream  temperatures  and 
temperature  changes,  and  monitored  the  particle  cloud  concentration. 

The  potentiometer  also  served  as  a  calibration  standard  for  the  other 
recorders . 


Energy  Balance  Procedures 

Energy  Balance  at  the  Entrance  of  the  Heat  Exchanger 

The  energy  absorbed  by  the  aerosols  was  determined  from  one  of 
two  energy  balance  calculations.  One  balance  was  made  between  the 
reference  state  and  the  thermal  conditions  at  the  entrance  of  the  heat 
exchanger  and  the  other  between  the  reference  state  and  the  exit  condi¬ 
tions.  They  were  made  by  taking  the  reference  temperature  as  that  of 
the  stream  at  the  exit  of  the  aerosol  generator  which  was,  effectively. 
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the  same  as  the  temperature  of  the  aerosol  at  the  entrance  of  the  furnace. 
This  hase  temperature  was  assumed  constant.  Actually,  variations  in  it 
never  exceeded  one  or  two  degrees  over  a  period  of  several  hours  and 
were  not  more  than  a  few  tenths  of  a  degree  over  the  duration  of  a  test. 

With  the  reference  temperature,  the  temperature  at  the  entrance 
of  the  heat  exchanger,  and  a  knowledge  of  the  properties  of  the  aerosol 
particles,  energy  balances  over  this  part  of  the  heat  transfer  system 
could  be  established.  The  difference  of  the  enthalpy  gain  per  unit  time 
of  the  air  in  the  stream  with  particles  and  the  gain  without  particles 
is  described  by 


“lA  =  '3.1) 

where  the  subscripts  (A,0,l)  and  the  superscripts  (O^l)  indicate  the 
following : 

A  =  a  property  of  the  air, 

0  =  the  property  evaluated  at  the  furnace  inlet, 

1  =  the  property  evaluated  at  the  heat  exchanger 
inlet , 

0  =  evaluation  when  no  particles  flowing, 

1  =  evaluation  when  particles  are  flowing. 

With  the  assumption  of  a  constant  base  temper^-ture,  equation  (3.l)  sim¬ 
plifies  to; 


^^lA  \a~\a 


(3.2) 


which  in  terms  of  temperature  is 
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where  C  is  the  average  heat  capacity  of  the  air  stream  hetween  the 

Jt\ 

temperatures  under  consideration^  and  is  the  mass  flow  rate  of  air. 

The  enthalpy  gain  per  unit  time  of  the  particles  in  the  streajn 
is  given  hy  the  relation 


AH,  =  hJ- 

Ip  Ip  Op 


(3.4) 


which  may  also  be  written 

“ip  - 

where  the  subscript  p  indicates  a  property  of  the  particles,  C  is  the 
average  heat  capacity  of  the  particle  material  between  the  temperatures 
under  consideration,  and  is  the  powder  flow  rate. 

The  total  enthalpy  gain  per  unit  time  due  to  the  presence  of 
particles  in  the  system  is,  finally,  the  srun 


AH,  .  , 
total 


+  AH, 
IP 


(3.6) 


Energy  Balance  at  the  Exit  of  the  Heat  Exchanger 

A  thermocouple  at  the  upper  end  of  the  heat  exchanger  sensed  the 
temperature  of  the  aerosol  stream  there,  and  a  difference  thermocouple 
indicated  the  temperature  rise  of  the  cooling  air  as  it  passed  through 
the  heat  exchanger.  From  these  temperatures  and  the  flow  rates  of  the 
streams,  a  second  energy  balance  was  established  that  included  the  heat 
exchanger.  The  total  balance  involved  the  sum  of  the  heat  effects  in 
the  air  portion  of  the  aerosol,  the  particles,  and  the  heat  exchanger 


coolant . 


For  air,  the  difference  of  the  enthalpy  gain  per  unit  time  with 
the  particles  and  the  gain  without  particles  is  given  by 

^3A  “  ~  (3*7) 

where  the  subscript  3  indicates  the  property  at  the  conditions  of  the 
heat  exchanger  exit. 

Since  the  base  temperature  was  constant,  equation  (3.7)  reduces 
to 

“3A  ‘  <3.3) 


or 


"•SA  ■  <“')A<‘k‘°3A> 


The  enthalpy  gain  of  the  particles  per  unit  time  is  then 


AH_  = 

3p  3p  Op 


(3.9) 


(3.10) 


which  can  be  expressed 


“3P  -  <*)p<4-‘^> 


(3.11) 


Finally,  the  difference  in  enthalpy  gain  per  unit  time  for  the  cooling 
air  is 


AH^  =  (WC)^  [(AtJ-At°)] 


(3.12) 


or, 


AH^  =  (WC)^  .  A(At^) 


(3.13) 
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where ^ 


A(At^)  =  [At^-At°]  (3.14) 

The  notation,  A(At^),  indicates  the  difference  in  temperature  rise  of 
the  heat  exchanger  coolant  with  and  without  particles  in  the  aerosol 
stream.  The  total  enthalpy  gain  per  unit  time  due  to  the  flow  of 
particles  is  thus  the  sum 

AH ,  ,  =  AH-„  +  AH-  +  AH  (3-15) 

total  3A  3p  c 

Materials  for  Study 

Selection 

Initial  measurements  were  made  with  zinc  powder  aerosols  since 
this  material  was  used  in  developing  the  apparatus  and  in  establishing 
experimental  procedures.  Being  a  metal,  it  has  a  low  absorptivity  and 
afforded  an  operational  test  of  the  apparatus  with  that  type  of  powders. 
Also,  the  zinc  was  well  chosen  for  testing  the  bag  filter  since  it  had 
an  average  particle  diameter  of  about  6  microns. 

Most  of  the  later  experimental  work  was  conducted,  however,  with 
ferrous  sulfide  and  copper  oxide  aerosols.  These  materials  were  selec¬ 
ted  primarily  because  of  their  noncorrosive  nature  and  their  availability 
in  lump  or  granular  form  which  permitted  them  to  be  subdivided  into 
various  size  fractions. 

Preparation 

The  zinc  powder  was  a  paint  pigment  obtained  from  the  Decatur 
Chemical  Company,  Decatur,  Georgia.  It  was  used  as  received.  Prolonged 
circulation  of  it  through  the  experimental  apparatus  during  the  period 
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of  development  eroded  the  individual  particles  until  they  became  almost 
spherical.  The  particle  diameters  then  ranged  from  one  micron  to  a 
maximum  of  10  microns  with  an  average  of  very  nearly  6.5  microns.  The 
actual  particle  distribution  curve  is  given  in  Figure  43  of  Appendix  A. 

The  ferrous  sulfide  was  the  technical  grade  product  of  the  Baker 
Chemical  Company,  Phillipsburg,  New  Jersey.  Preparation  of  it  consisted 
of  first  reducing  the  lumps  to  minus  ikk  U.S.  Standard  Sieve  Size  with  a 
mortar  and  pestle.  The  product  was  then  ball-milled  for  one  hour  inter¬ 
vals  and  screened  until  all  of  it  passed  an  88  micron  screen.  There¬ 
after,  standard  sieving  was  continued  to  yield  88  to  53 >  53  to  44,  and 
less  than  44  micron  diameter  fractions.  Special  screens  (Buckbee 
Mears  Co.,  St.  Paul,  Minn.)  were  then  employed  to  separate  the  less 
than  44  micron  material  into  44  to  30,  30  to  20,  and  less  than  20  micron 
diameter  fractions. 

The  cupric  oxide  was  also  reagent  grade  material  and  obtained 
from  the  Baker  Chemical  Company.  This  high  purity  product  was  employed 
since  only  it  could  be  located  with  particles  larger  than  20  microns. 
Particle  size  fractions  ranging  from  the  88  micron  screen  size  to  20 
microns  were  prepared  by  the  same  procedure  as  was  the  ferrous  sulfide. 
Particle  Properties  and  Size  Measurements 

After  preparation  of  the  different  fractions,  precise  measure¬ 
ments  of  particle  diameters,  distribution  of  diameters,  and  surface 
area  were  made,  since  these  parameters  were  thought  to  be  significant 
in  establishing  the  absorptive  properties  of  particle  clouds  formed 
from  them.  Attempts  were  also  made  to  establish  the  size  of  cracks 
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and  fissures  within  the  particles.  These  tests  involved  microscopic 
determinations  of  thousands  of  particle  diameters,  statistical  treat¬ 
ment  of  the  data,  and  low  temperature  nitrogen  absorption  measurements 

40 

in  accordance,  with  standard  procedures.  Effective  particle  surface 
areas  were  calculated  from  the  size  distributions,  for  it  was  this 
quantity  that  determined  the  specific  absorptive  capability  of  the 
powder  (radiation  properties  excepted).  It  was  anticipated  that  the 
nitrogen  absorption  surface  areas,  when  allowance  was  made  for  the 

4l  42 

porous  structure  of  the  surface,  ^  might  be  employed  to  yield  an 

2 

effective  specific  surface  area  (ft  /lb).  The  attempts  might  have 
been  successful  with  some  powders,  particularly  with  very  small,  hard- 
to-disperse  particles  where  the  nitrogen  absorption  measurements  are 
more  reliably  made  than  are  microscopic  determinations.  Unfortunately, 
adsorption  determinations  could  not  be  obtained  with  adequate  preci¬ 
sion.  There  still  remains  the  possibility  that  krypton  adsorption 

43 

might  make  the  method  reliable.  More  research  is  required,  however, 

* 

before  pore  structure  can  be  determined  with  krypton. 

Physical  measurements  were  made  both  before  and  after  several 
of  the  powder  fractions  were  used  in  the  heat  transfer  experiments  to 
determine  if  any  changes  in  properties  had  resulted.  As  expected, 
significant  size  and  shape  differences  were  found,  particularly  in  the 
first  experiments  that  required  long-time  operations  to  produce  useful 
results.  Prolonged  particle  attrition  produced  broadened  size 

^  — 

This  is  being  pursued  in  another  thesis  investigation  now  being 
conducted  by  Korstiaan  Van  Wijngaarden  of  this  laboratory. 
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distributions  and  resulted  in  particles,  in  some  cases,  well  below  the 
intended  size.  The  powders,  appeared  to  change  most  rapidly  during 
the  first  few  minutes  of  use*  Thereafter  changes  were  much  slower, 
and  size  distributions  determined  immediately  after  use  were  employed 
for  the  correlations.  These  distributions  are  given  in  Appendix  A, 

One  of  the  most  unexpected  observations  is  presented  in  Table  1. 
There  it  may  be  seen  that  the  powders’  specific  surface  areas  were 
reduced  by  use.  The  fact  at  first  seemed  paradoxical  since  the  parti¬ 
cles  undoubtedly  underwent  size  reduction.  The  effect,  however,  resul¬ 
ted  from  the  capture  in  the  bag  filter  of  submicron  particles  which 
initially  were  present  but  were  attached  to  the  large  particles.  These 
were  not  included  in  the  original  particle  size  distribution  measure¬ 
ments.  Figure  25  is  a  series  of  photographs  of  a  powder  before  use, 
the  portion  retained  in  the  bag,  and  a  fraction  after  use. 

Calculated  Particle  Cross  Sectional  Areas 

The  particle  area  actually  participating  in  the  radiation  inter¬ 
change  was  evaluated  from  an  experimentally  determined  quantity,  S^, 
the  cross  sectional  area  per  unit  volume  of  solids  (ft  was 

constant  for  each  powder  studied  and  was  a  function  of  particle  geo¬ 
metry  and  particle  size  distribution.  Its  value  was  calculated  as 
follows : 


For  a  spherical  particle  of  diameter. 


d  ,  the  volume  is 

p 


(3.16) 


and  the  cross  sectional  area  is 
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Table  1.  Surface 


53-88 

44-53 

30-44 


jas  of  Powdered  Materials  as 
by  Nitrogen  Absorption 


Specific  Surface  Area 

Before  Use 

After  Use 

(m^/gm) 

(m^/gm) 

0.37 

0.13 

0,35 

0.l4 

0.45 

0.23 

(a)  Original  Material 


(b)  Material  Retained  in  Bag  Filter 


(c)  Circulated  Material 


Photomicrographs  of  Ferrous  Sulfide  Powder 
Showing  the  Preferential  Capture  of  the  Very  Small 
Particles  in  the  Bag  Filter. 


A  =  nd^ 
c  p 


(3.17) 


The  ratio  of  to  V  gives  which  in  this  case  is  6/d^.  For  a 
powder  having  a  range  of  particle  sizes ^  the  contributions  from 
particles  within  a  differential  diameter  increment  may  be  written 

■  <s>(diy)(''c)nd 

and 

■i'"  (3-W) 

which,  when  integrated  give  the  total  area  and  total  particle  voluime  of 
N  particles,  respectively 

'  f  '>(%>  (3.20) 

P 

and 

"■f  (s)(d^)(=r>^("p)  (3-^^) 

1  ^ 

Again  is  the  ratio  A^/V. 

Under  certain  conditions  the  indicated  integrations  have  clas¬ 
sical  solutions.  If,  for  instance,  the  particle  size  distributions 

44 

were  log-normally  distributed  the  two  results  would  be 

Aq  =  -^  Ag  =  .^  exp  (2  In  Xg  +  2  In^  a^)  (3-22) 

V  =  exp  (3  In  X  +  4.5  In^  a  )  (3.23) 

D  D 
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Where  x  is  the  geometric  mean  diameter  and  a  is  the  geometric  standard 
deviation.  Many  powders  in  nature  can  be  described  accurately  by  this 
distribution  and  for  such  cases  the  final  result  for  is  immediately 
available . 

The  powders  employed  in  this  investigation,  unfortunately,  did 
not  follow  the  log  normal  distribution  law,  since  this  type  of  distribu¬ 
tion  was  destroyed  by  the  sizing  processes.  Log  normal  particle  size 
distributions  are,  in  general,  produced  by  natural  comminution.  The 
actual  particle  size  distributions,  in  some  cases,  were  bimodal. 

Attrition  in  the  heat  transfer  apparatus,  if  continued,  would  undoubtedly 
eventually  have  produced  a  log  normal  distribution  but  much  less  than 
this  length  of  time  was  allowed. 

A  numerical  procedure  was  employed  to  evaluate  the  surface-to- 
particle  volume  ratio  since  no  mathematical  representation  for  describ¬ 
ing  the  distribution  was  available.  Equations  (3.20)  and  (3-21),  were 
integrated  numerically  by  taking  diameter  increments  of  about  one  to 
two  microns  and  computing  the  required  areas  and  volumes  for  the  average 
diameter  of  the  increment.  The  results  were  summed,  and  area- to- volume 
ratio  were  formed  to  give  the  final  answer.  These  calculations  were 
performed  with  the  Burroughs  220  digital  computer. 

Experimental  Data 
With  Water-Cooled  Heat  Exchanger 

The  heat  transferred  to  a  given  dispersion  was  determined  by 
energy  balances  as  already  described.  For  the  initial  experiments, 
water  was  the  coolant  for  the  heat  exchanger  and  the  energy  balance 
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given  ty  equation  (3.15)  was  employed  to  calculate  the  heat  absorbed  by 
the  aerosol.  Data  were  collected  in  this  manner  for  the  zinc  powder 
and  the  53-88  micron  diameter  fraction  of  ferrous  sulfide.  T^ese 
experimental  and  calculated  results  are  given  in  Tables  3  and  4  of 
Appendix  B. 

The  tabulated  results  for  these  experiments  include  the  thermal 
data  for  the  streams  involved;  the  air  and  powder  flow  rates;  the 
indicated  furnace  temperature;  and,  finally,  the  calculated  enthalpy 
changes  produced  by  the  addition  of  particles  to  the  system.  The 
temperatures  for  the  aerosol  and  the  heat  exchanger  coolant  are  pre¬ 
sented  directly  as  indicated  by  thermocouples  located  as  shown  on 
Figure  17  (page  65).  The  difference  quantity,  A(At),  however,  is  a 
combination  of  two  measurements,  one  the  temperature  rise  of  the  heat 
exchanger  coolant  when  only  air  was  flowing  in  the  system,  and  the 
other  when  both  air  and  powder  were  flowing.  The  temperature  rise  of 
the  heat  exchanger  coolant  was  a  temperature  difference,  hence  the 
notation  (At).  The  difference  between  these  two  (At)  quantities 
represents  the  change  in  temperature  rise  of  the  heat  exchanger 
coolant  when  particles  were  added  to  the  system.  The  second  differ¬ 
ence  was  represented  by  A(At).  This  quantity  in  some  cases  was  posi¬ 
tive  and  in  others  negative,  depending  on  whether  more  or  less  heat 
was  absorbed  by  the  heat  exchanger  when  particles  were  added  to  the 
system. 

Making  an  energy  balance  for  these  experiments  required,  as 
shown  by  equation  (3.15)?  the  enthalpy  changes  incurred  by  the  aerosol 
components  as  well  as  the  heat  exchanger  coolant  when  particles  were 
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flowing.  The  latter  could  not  be  calculated  directly  but  was  estimated 
to  be  comparatively  small.  Tables  3  and  4,  therefore,  present  only  the 
calculations  for  the  aerosol  particles  and  the  carrier  air. 

As  shown  by  the  tabulated  data,  the  A(At)  for  the  heat  exchanger 
coolant  was  approximately  ±  0.5°  F.  With  a  water  flow  rate  of  0.095  CFM 
this  means  a  (AHc)  of  3  Btu/min.  The  result,  however,  was  unreliable 
since  the  possible  error  in  determining  A (At)  was  on  the  order  of  100 
per  cent.  Other  measurements,  therefore,  were  required  to  assess  this 
quantity  accurately.  A  thermocouple  was  installed  in  the  aerosol  stream 
near  the  entrance  of  the  heat  exchanger,  as  previously  described,  and 
was  employed  as  a  supplementary  measurement  in  analyzing  the  operations 
of  the  heat  exchanger.  Since,  with  the  new  measurement,  the  entrance 
and  exit  conditions  of  both  streams  were  known,  log-mean  temperature 
differences  between  the  coolant  and  the  aerosol  system  could  be  calcu¬ 
lated.  These  differences,  when  powder  was  flowing,  were  found  to  range 
from  5  to  15  per  cent  less  than  these  differences  when  only  air  was  flow¬ 
ing.  It  was  noted  further  that  the  overall  heat  transfer  coefficient 

45 

was  correspondingly  increased  for  the  powder  concentrations  studied. 

These  conditions,  then,  meant  that  the  heat  transferred  (UA  At)  to  the 
coolant  in  each  case  was  essentially  the  same,  and  it  followed  further 
that  the  A(At)  must  have  been  very  nearly  zero..  Thus,  it  was  assumed 
that  the  heat  exchanger  absorbed  approximately  the  same  amount  of  heat 
with  and  without  powder  and  the  enthalpy  change  AHc  for  the  coolant  was 
disregarded.  These  results  are  shown  in  Figure  26  for  the  zinc  and  in 
Figure  27  for  the  53  to  88  micron  ferrous  sulfide  particle  aerosols. 
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Subsequent  data  with  other  ferrous  sulfide  powders,  also  presented  in 
the  latter  figure,  gave  support  to  this  conclusion. 

With  Air-Cooled  Heat  Exchanger 

For  optimum  operation  the  heat  exchanger  should  be  near  room 
temperature  to  reduce  to  a  minimum  extraneous  heat  effects  associated 
with  the  surroundings.  For  this  purpose  water  is  a  better  coolant  due 
to  its  greater  heat  capacity.  Its  use,  however,  was  actually  a  dis¬ 
advantage  in  evaluating  the  amount  of  energy  absorbed  by  the  water  with 
acceptable  precision.  When  the  water  rate  was  sufficiently  high  to 
provide  stable  operation,  the  temperature  differences  involved  were  too 
small  to  measure  accurately,  and  when  the  water  rate  was  reduced  to  give 
measurable  temperature  differences,  the  heat  exchanger  was  not  thermally 
stable  and  too  much  time  was  required  for  equilibration.  Air,  therefore, 
was  employed  as  the  heat  exchanger  coolant  in  an  effort  to  achieve  better 
measurements  and  more  stable  operation.  Air,  having  a  much  smaller  heat 
capacity  than  water,  could  have  its  flow  rate  regulated  to  give  measur¬ 
able  temperature  differences  as  well  as  rapid  response  to  transient 
heat  transfer  effects.  Data  were  collected  with  this  modification  for 
the  44  to  53?  30  to  44  and  the  20  to  30  micron  ferrous  sulfide  powders. 
The  results  of  these  experiments  are  presented  in  Tables  5  through  8  and 
on  the  preceding  figure.  The  absorption  determinations  were  reprodu¬ 
cible  within  ±  0.5  Btu/min.  and  showed  a  considerable  reduction  in  the 
uncertainty  of  the  heat  removed  by  the  heat  exchanger  coolant  as  com¬ 
pared  with  the  results  using  water.  Air  was  used  in  all  subsequent 
e}<periments.  Also,  starting  with  the  20  to  30  micron  ferrous  sulfide 
powder  (Table  8),  data  were  recorded  so  that  heat  balances  could  be 
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made  over  the  complete  heat  exchanger  [equation  (3*15)]  well  as  on 
the  aerosol  system  alone  [equation  (3.6)],  except  for  one  e:Kperiment 
with  cupric  oxide  (Table  12), 

Final  measurements  were  made  on  aerosols  of  cupric  oxide  powders. 
Four  size  distributions  were  prepared  as  previously  described.  Again 
samples  designated  by  screen  sizes  of  53  to  88,  44  to  53?  30  to  44  and 
20  to  30  microns  were  selected.  Tables  9  through  12  present  the  data 
obtained  while  Figure  28  gives  the  results  graphically. 

Radiant  Heat  Flux  Determinations 

The  radiant  heat  flux  penetrating  the  quartz  enclosures  and 
traversing  the  particle  cloud  was  basic  to  the  theoretical  heat  transfer 
calculations  of  this  study.  To  determine  this  quantity,  the  aerosol 
conduit  was  coated  with  a  highly  absorbing  material  --  carbon  black  -- 
and  experiments  were  conducted  in  which  water  was  forced  to  flow  upward 
through  the  inner  quartz  conduit  to  remove  the  heat  absorbed  thereon. 

The  contribution  of  the  radiation  had  to  be  separated  from  the  total 
result,  however,  because  the  aerosol  conduit  received  heat  by  conduc¬ 
tion  and  convection  as  well  as  radiation.  This  separation  was  accom¬ 
plished  by  performing  additional  experiments  with  the  inner  wall  of  the 
aerosol  conduit  being  made  as  reflective  as  possible  so  that  radiation 
effects  would  be  minimized.  A  gold  coating  was  employed  for  this  pur¬ 
pose.  The  difference,  then,  between  the  results,  when  adjusted  for  the 
actual  absorptivities  of  the  two  coatings,  was  the  heat  transferred  to 
the  aerosol  conduit  by  radiation. 

The  radiant  heat  flux  determinations  were  made  with  water  as  the 
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coolant,  rather  than  air,  to  permit  better  control  of  the  temperature  of 
the  system.  When  air  was  the  coolant  and  the  conduit  wall  was  made 
highly  absorptive,  a  sufficiently  high  flow  rate  of  air  could  not  be 
accommodated  by  the  rotameter  metering  the  air  to  avoid  thermal  damage 
to  the  coating.  Further,  it  was  desired  that  the  conductive  and  the 
convective  heat  transfer  to  the  system  be  as  nearly  identical  as  pos¬ 
sible  with  both  the  reflective  and  the  absorptive  coatings.  The  en¬ 
trance  and  exit  temperatures  of  the  coolant  had  to  be  maintained  suf¬ 
ficiently  close  to  avoid  large  discrepancies,  since  both  these  processes 
were  a  function  of  the  temperature  distributions  throughout  the  complete 
system.  With  water  as  the  coolant,  flow  rates  could  be  selected  to 
keep  the  temperature  rises  within  about  5*^  F  of  each  other.  For  a  pair 
of  experiments,  one  with  each  type  of  coating,  the  flow  rates  had  to  be 
identical  since  the  extent  of  the  convective  and  conductive  heat  trans¬ 
fer  was  a  function  of  the  coolant  flow  rate. 

The  conduit  wall  was  made  highly  absorbing  by  coating  it  with 
lamp  black  in  a  small  amount  of  epoxy  resin  (Resinweld  Plastic  Alloy, 

H.  B.  Fullter  Co.,  St.  Paul,  Minn.).  The  absorptivity  of  the  surface 
was  taken  as  The  reflective  surface  was  a  thin  gold  film  baked 

onto  the  inner  conduit  wall.  A  fluid  suspension  of  gold  particles 
supplied  by  B.  F.  Drakenfeld  and  Co.  Inc.,  Washington,  Penn,  was  used. 
Baking  a  deposit  of  this  suspension  to  a  temperature  (l060^  F)  near 
the  melting  point  of  gold  volatilized  the  suspending  fluids  and  formed 
a  smooth  reflective  metallic  surface.  The  absorptivity  was  assumed  to 

4? 

be  0.03. 

Experimental  and  calculated  heat  transfer  data  are  presented 
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graphically  in  Figures  29,  30,  and  31.  The  first  two  figures  give  the 
total  heat  absorbed  by  the  blackened  and  the  reflective  walls,  re¬ 
spectively,  as  a  function  of  the  indicated  furnace  temperature  and 
three  water  flow  rates.  The  latter  figure  presents  the  furnace  heat 
flux  as  a  function  of  temperature  after  corrections  were  made  for  the 
small  absorptivity  deviation  in  the  emissive  properties  of  the  wall 
coatings,  i.e.,  0.95  instead  of  1.00  for  the  blackened  wall  and  0.03 
instead  of  0.00  for  the  reflective  coating.  At  the  indicated  operating 
temperature  of  the  furnace  (2280  F),  the  heat  flux  was  17,600  Btu/hr 


CHAPTER  IV 


RESULTS  MB  DISCUSSION 


Theoretical  Developments 
Generalized  Calculations 

The  theoretical  considerations  of  this  study  led  to  the  develop¬ 
ment  of  an  equation  for  evaluating  the  direct  radiant  heat  transfer 


from  a  hlack  cylindrical  enclosure  to  a  suspension  of  "black  particles 


contained  within.  The  final  result,  as  given  previously  by  equations 

(2.57)  and  (2.58), is 


P,  =2ni 
b 


^nw 


CS  R 

V 


(Jl-vJj.) 


rdr 


(4.1) 


and  with  the  integral  evaluated 


Q,  =2)ti  CS  R(f(+F') 
b  Tjnw  V  '  1  E'^ 


(4.2) 


These  relationships,  however,  are  dimensional;  they  may  be  transformed 
into  a  more  useful  form  by  rendering  them  dimensionless.  Thus, 


and 


=  2kR 


(Jl+J^)  (j^)d(j^) 


(4.3) 


- - ^  =  2kR(F  +F_) 


(4.4) 


where  F^,  F^,  and  k  are  F^/R^,  F^R^,  and  CS^  respectively.  The  rate 
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of  heat  absorption  (Q^)  ^'3.  the  product  (i^^^)(TTR  )  have  the  units  of 
Btu  per  unit  time;  their  ratio  is  therefore  non-dimensional.  The 
quantity,  KR,  is  the  absorption  coefficient  (units  =  reciprocal  length) 
times  the  radius  of  the  system  enclosure  and  is  therefore  a  characteris¬ 
tic  quantity  independent  of  the  length  unit.  From  equations  (2.19), 
(2.43),  (2.57)  and  (2.58)  it  may  be  seen  that  F^  is  a  function  of  L/R 
only  and  F^,  is  a  function  of  L/R  and  kR. 

A  computer  program  for  the  Burroughs  5000  digital  computer  of 
the  institution's  Rich  Electronic  Con^uter  Center  was  developed  for 
evaluating  the  terms  of  the  heat  transfer  equation  and  a  general  study 
of  it  was  made.  Sufficient  results  were  calculated  to  construct  a 
generalized  plot  describing  the  radiation  process.  The  results  are 
given  in  Figure  32  where  the  quantity  plotted  versus 

L/R  with  kR  as  a  parameter.  The  ordinate  represents  the  rate  of  heat 
absorbed  per  unit  radiation  intensity  per  unit  cross  sectional  area  and 
is  solely  determined  by  the  geometrical  term  L/R  and  the  absorption 
parameter  kR.  For  a  constant  L/R,  the  absorption  ordinate  increases 
with  the  absolute  value  of  kR  and  finally  when  kR  is  infinite  the 
absorption  is  a  maximum  given  by  the  uppermost  line.  This  limiting 
curve  represents  the  total  emissive  power  of  the  system  enclosure,  and 
the  equation  for  it  is 


%mit _ 


2n(i) 


(“i.S) 
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Being  that  the  total  emitted  energy  is  represented,  the  plot  has  the 
desirable  feature  that  the  absorption  efficiency  is  immediately 


_ Heat  Absorbed _ 
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Figure  32.  Generalized  Correlation  for  the  Absorption  of  Radiant 
Energy  by  an  Evenly  Dispersed  Particle  Cloud  within  a 
Unifonnly  Radiating  Cylindrical  Enclosure. 
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available  by  inspection;  it  is  only  necessary  to  divide  the  ordinate  for 
a  given  system  by  the  total  emission  represented  by  the  uppermost  curve. 
The  efficiency,  then  is 


Y  = 


Q 


emit 


^emit 


(4.6) 


Each  parametric  curve  has  three  distinct  regions  classifiable  by  the 
L/R  ratio.  For  very  small  or  very  large  values  of  the  abiscissa,  each 
curve  is  very  nearly  a  straight  line.  Between  these  regions  each 
curve  has  a  variable  slope  such  that  the  functions  are  smoothly  contin¬ 
uous  between  the  two  characteristically  linear  regions.  The  linear 
extensions  of  the  curves  make  them  suitable  for  indefinite  extrapolation. 
Construction  of  Additional  Parametric  Curves 

While  the  generalized  graph  gives  accurately  the  effect  of  the 
geometric  variables,  the  particle  cloud  absorption  properties  were 
limited  to  selected  parametric  curves  for  kR,  and  interpolation  between 
them  is  uncertain.  The  data,  however,  may  be  replotted  to  yield  results 
for  any  value  of  kR,  if  L/R  is  made  the  parameter.  The  graph,  for  the 
most  part,  is  very  nearly  linear,  except  near  total  absorption,  and 
results  for  a  particular  kR  can  be  read  very  accurately.  If  only  one 
system  is  being  studied,  the  plot  is  made  for  the  required  L/R.  If, 
however,  a  completely  new  parametric  curve  is  needed  the  process  is 
repeated  for  several  L/R  values  and  as  many  points  as  desired  can  be 
plotted  on  the  generalized  graph,  and  the  new  curve  constructed.  Five 
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charts  of  this  nature  are  given  in  Figures  33  through  37  for  several 
widely  spaced  L/R  values.  The  construction  of  additional  curves  by 
this  method  is  very  accurate  and  results  from  the  procedure  were 
checked  with  the  computer  and  almost  exact  agreement  was  found.  For 
points  near  total  absorption,  however,  calculated  data  were  not 
obtained  because  of  arithmetic  or  exponential  overflows  in  the  com¬ 
puter  operations,  unsuitable  accuracy  in  partial  results,  and  \msat- 
isfactory  convergence.  Extrapolation  of  these  data  to  kR  values 
larger  than  those  given  may  be  attempted  but  the  results  will  be 
outside  the  range  of  computations. 

Use  of  the  Generalized  Chart 

The  generalized  graph  was  developed  specifically  for  the 
evaluation  of  the  absorption  of  radiant  energy  by  clouds  of  uniformly 
dispersed  black  particles  while  within  a  black,  uniformly  radiating 
cylinder.  The  result  for  a  particular  problem  is  obtained  by  forming 
the  length-to-radius  ratio  and  the  absorption  parameter,  kR;  deter¬ 
mining  the  corresponding  ordinate  for  these  values  and  multiplying 
this  number  by  the  radiation  intensity  and  the  cross  sectional  area 
of  the  enclosure.  The  result  is  the  heat  absorbed  by  direct  emission 
from  the  radiating  wall;  the  value  divided  by  the  total  emission  from 
the  enclosure  gives  the  particle  cloud  absorption  efficiency.  This 
is  the  primary  use  of  the  generalized  graph. 

If  the  total  length  of  an  enclosure  is  not  heated  and  unheated 
sections  extend  beyond  the  heated  region,  there  is  additional  absorp¬ 
tion  of  radiation  by  particles  contained  in  the  unheated  extensions. 
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Figure  3h,  Generalized  Results  for  a  Length-to-Radius  Ratio  of  0.1. 
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Figure  35-  Generalized  Results  for  a  Length-to-Radius  Ratio  of  1. 
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Figure  36.  Generalized  Results  for  a  Length-to-Radius  Ratio  of  10.0. 


o 


Figure  37 •  Generalized  Results  for  a  Length-to-Radius  Ratio  of  48.906. 


This  absorption  may  be  significant  in  some  systems,  especially  when  the 

length-to-radius  ratio  is  small.  The  increase  in  heat  transfer  is 

calculable  from  the  basic  radiation  equation  by  changing  the  limits  of 

integration  and  developing  a  calculation  procedure  similar  to  the  one 

presented.  The  same  result,  however,  is  available  from  the  generalized 

chart,  but  the  following  developments  are  required. 

The  absorption  process  in  an  unheated  length  can  be  described 

by  first  analyzing  how  absorption  within  a  heated  region  is  changed 

upon  adding  other  heated  increments.  Consider  two  equal  lengths  of 

enclosure.  If  they  are  adjacent,  the  heat  absorbed  may  be  designated 

as  E  ;  if  they  are  separate  the  total  for  both  is  2E  .  The  absorption 
2x’  X-  X 

is  greater  when  the  two  lengths  are  adjacent  because  of  mutual  absorp¬ 
tion  from  each  other.  It  follows  then  that 

E^  -  2E  =  2AE^  (i+.Y) 

2x  X  X 

where  2AE^  is  the  increase  in  absorption  when  the  lengths  are  placed 

X 

together.  The  subscript  on  this  quantity  indicates  the  absorbing  length 
and  the  superscript  indicates  the  emitter.  The  factor  2  appears  because 
of  the  reciprocal  emission  and  absorption  between  the  two  lengths. 

The  next  step  is  to  consider  three  increments.  The  relationship 
for  this  case  is 

E.,  -  3E  =  2AE^’^  +  2AE^  (4.8) 

where  again  the  delta-quantities  describe  the  effect  of  placing  the 
three  lengths  together.  They  are  derived  by  considering  the  radiation 
received  by  each  length,  taken  individually,  that  is  emitted  from  all 
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lengths  on  each  side  of  it^  taken  separately.  By  way  of  explanation, 

define  the  increments  consecutively  as  A,  B,  and  C.  If  A  and  B  are 

taken  together  radiating  to  C,  the  result  is  AE^  (here  2x  is  the 

emitting  length  and  x  the  absorber).  The  lengths  B  and  C  radiate  to 

2x  X 

A,  hence  a  second  AE  .  The  total  then  is  2AE  .  Then  A  can  radiate 

’  X  X 

to  B  and  C  to  B;  these  combinations  give  2AE  . 

If  four  increments  (A,  B,  C,  and  D)  are  taken,  a  similar  equa¬ 
tion  with  an  additional  term  is  obtained.  Again  the  radiation  from 
the  entire  neighborhood  to  each  element  is  considered;  a  different 
approach  is  clearer,  however,  when  higher  combinations  are  involved. 
The  notation  may  be  employed: 


ABCD  -*  A  =  (a 

ABCD  -*  B  =  (B 

ABCD  -*  C  =  (C 

ABCD  D  =  (D 

and 

(ABCD  -  A)  +  (ABCD  ^  B)  + 


-*  a)  + 

(BCD 

-  A) 

(4.9) 

-*  b)  + 

(acd 

-  B) 

(4.10) 

c)  + 

(abd 

-c) 

(4.11) 

-*  D  )  + 

(ABC 

-  D) 

(4.12) 

(ABCD  - 

-  C)  -t 

■  (ABCD  -  D)=Ei^^ 

(4.13) 

where  the  arrow  means  "radiates  to."  If  these  results  are  added  to- 
gether  the  total  radiation  recieved  by  the  elements  in  When  the 

elements  are  separate  the  radiation  interchange  may  be  defined  as 


A  -*  A  (4.l4) 

B  -  B  (4.15) 


C  ->  C  (4.l6) 

D  D  (4.17) 
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and 


(A  -»  A)  +  (B  -  B)  +  (C  -  C)  +  (D  -  D)  =  i+E  (4.l8) 

Here  the  total  radiation  absorption  is  4E  .  Subtraction  of  the  latter 
set  of  relations  from  the  former  gives 


(ABCD  -  A) 

_  (a  -  a)  =  BCD  -  a 

(4.19) 

(abcd  b) 

-  (B  B)  =  ACD  B 

(4.20) 

(abcd  -  c) 

-  (C  -*  C)  =  ABD  -  C 

(4.21) 

(abcd  -*  D) 

-  (D  -*  D)  =  ABC  -*  D 

(4.22) 

Following  the  development  already  given  for  equation  (4,8)  it  may  be 
shown  that 


BCD 

-.A.< 

(4.23) 

ABC 

D  =  AE^^ 

X 

(4.24) 

ACD 

B  =  AE^^  AE^ 

X  X 

(4.25) 

ABD 

-  C  =  AE^""  +  AE^ 

(4.26) 

Equations  (4.9)  through  (4,26)  may  be  shown  to  combine  and  give 

E,  -  4e  =  2AE^^  +  2AE^^  +  2AE^  (4.2?) 

4x  X  X  X  X 

If  more  increments  are  added,  a  series  is  established.  In 
general  then,  it  follows  that 

E/  -  (n+l)E  =  2AE’^^  +  2AE^’^"^^^  +  ...  +  2AE^  (4.38) 

(n+l)x  ^  ^  X  X  X  X 
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and 


E  -  nE  =  +  ...  +  2AE^  (4.29) 

nx  X  X  X  X 

These  equations  may  be  rearranged  and  solved  to  give 

AE™  =  ^[  E,  -  (n+l)E  ]  -  [  AE^^”^^^  +  AE^^“^^^  +  ... 

X  2  (ii+l)x  ^  X  ■’  X  X 

+  AE^  ]  (4.30) 

X. 


and 


AE 


nx 

X 


•ICVDx-  -5 


nx 


nE  ) 


(4.31) 


Since 


(n+l)E^ 


nE  +  E  , 

X  X 


(4.32) 


AE 


nx 

X 


^(n+l)x 

2 


E 

nx 


E 

X 


(4.33) 


The  significance  of  AE^  is  that  it  arises  from  the  end  incre¬ 
ment  receiving  radiation  from  all  the  n  other  elements.  This  amount 
of  energy  is  received  regardless  of  whether  or  not  the  end  element  is 
emitting.  Hence,  if  nx  is  a  heated  length,  the  radiation  from  it  to 
X  can  be  evaluated  from  equation  (4.33)« 

When  a  length  of  enclosure  radiates  to  a  longer  length  it  is 
X  nx 

desirable  to  have  AE  rather  than  AE  .  The  result  is  directly 

nx  x 

available  from 


•n-y  -y- 

E,  ,  X  -  E  -  E  =  AE  +  AE 
(n+l)x  nx  XX  nx 


which  together  with  Equation  (4.33)  yields 


(4.34) 
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(4.35) 


AE^ 


nx 


-  E 

nx 

2 


E 

X 


It  follows  also  that 


AE 


X 

nx 


(4.36) 


The  derivations  are  restricted  to  n  being  an  integer.  It  may 
be  readily  shown,  however,  that  the  derivation  holds  for  any  real 
value  of  n. 

The  utility  of  the  extension  formulas  may  be  demonstrated  by 
specific  examples.  Take,  for  instance,  a  cylinder  with  an  L/R  ratio 
of  3*0  containing  a  particle  cloud  with  a  kR  value  of  1.0.  When  an 
unheated  extension  equal  to  one  half  the  heated  length  is  involved, 
equation  (4.33)  is  applicable  and  n  is  2.  The  heat  absorbed  by  parti¬ 
cles  while  they  are  in  the  extension  is 


(4.37) 


where  is  the  generalized  ordinate  for  L/R  of  1.3,  for  an  L/R 


X 


of  3.0,  and  E_  for  an  L/R  of  4,5,  each  being  established  on  the 
3x 

curve,  kR  =  1.  The  points  for  the  L/R  ratios  reflect  the  changes  in 
cylinder  length  being  considered,  since  the  radius  is  the  same  for 
all  cases.  If  there  are  extensions  on  both  ends  of  the  heated  length, 
calculations  are  made  for  each  separately.  When  both  extensions  are 
the  same,  the  result  for  one  length  is  doubled. 

In  the  preceding  example,  if  the  unheated  length  were  five  times 
that  of  the  heated  section  the  result  would  be 
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(4.38) 


—  E  —  E 

AT7^  -  OX  X 

where  again  x  is  the  L/R  ratio  of  the  basic  heated  length.  The  unheat- 
ed-to-heated  length  ratio  does  not  have  to  be  integral.  For  instance, 
it  may  be  5.5  and  the  additional  heat  absorption  would  be 

^  (‘*•39) 

X  nx 

For  each  example,  AE  ,  or  AE  ,  whichever  is  needed,  is  to  be  added  to 
E^  for  the  heated  length  to  obtain  the  total  result  for  the  system. 

In  performing  the  calculations  for  specific  cases  it  may  be  foimd 
that  the  total  absorption  (E  +  AE)  does  not  increase  smoothly  with  in¬ 
creasing  unheated  length.  This  condition  results  from  inaccuracies  in 
reading  the  generalized  graph.  To  obtain  improved  results,  it  may  be 
necessary,  in  particular  for  L/R  over  unity,  to  construct  a  larger,  more 
accurate  graph  by  the  use  of  Figures  33  through  37  as  already  described. 
Application  to  Other  Types  of  Radiation 

While  this  investigation  was  concerned  with  thermal  radiation, 
the  basic  equations  and  results  are  applicable  to  all  forms  of  electro¬ 
magnetic  and  field- emission  radiation.  One  conceivable  use  might  be 
the  evaluation  of  the  radiation  dosage  an  aerosol  receives  while  passing 
through  a  cylindrical  enclosure  with  nuclear  radiation  coming  into  its 
interior.  The  calculation  may  be  made  directly  if  a  radiation  cross- 
section  value  for  the  suspended  particles  is  known. 

Other  Important  Variables 

For  many  systems,  the  described  calculations  are  adequate.  Fre¬ 
quently,  however,  either  the  wall  enclosure  or  the  particle  surfaces. 
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or  ’both,  are  not  perfect  absorbers  and  adjustments  in  the  calculations 
may  be  required.  The  emissivity  of  the  wall  is  first  taken  into  account 
when  the  radiation  intensity  is  evaluated;  the  calculations  as  described, 
however,  are  unaffected  at  this  point.  If  radiation  crosses  the  enclo¬ 
sure  and  is  intercepted  by  the  wall  surface,  it  is  partly  absorbed  and 
partly  reflected.  The  reflected  energy  is  available  for  further  absorp¬ 
tion  by  the  particles  and  the  infinitely  continuing  process  adds  to  the 
heat  transfer  as  more  reflections  occur.  Only  the  first  few,  perhaps 
two  or  three,  need  be  considered  in  many  cases. 

If  the  energy  intercepted  by  the  wall  from  each  emission  is 
evaluated,  these  reflections  may  be  traced  and  their  contribution  to 
the  overall  result  calculated.  It  is  feasible  to  attempt  these  calcula¬ 
tions,  but  a  detailed  study  is  required  before  general  conclusions  can 
be  stated.  It  is  reasonable  to  expect,  however,  that  situations  with 
large  wall  emissivities  should  correspond  very  closely  to  the  ideal 
solution.  If,  for  example,  the  wall  emissivity  is  0.75  then  only  one- 
fourth  of  the  radiation  intercepted  by  the  walls  is  reflected.  If  the 
aerosol  is  highly  absorbing  very  little  energy  reaches  the  wall,  hence 
the  reflected  energy  may  be  neglected.  If  the  aerosol  is  weakly  absorb¬ 
ing,  a  large  portion  of  the  energy  may  be  intercepted  by  the  wall,  the 
extent  of  interception  being  dependent  on  the  wall-to-wall  view  factor. 
When  the  reflected  energy  traverses  the  aerosol,  part  of  it  Xs  absorbed 
therein,  and  in  terms  of  a  percentage,  this  part  is  approximately  the 
same  as  for  the  original  direct  emission.  Since  a  maximum  of  25  per 
cent  of  the  direct  emission  can  be  reflected,  the  additional  absorption 
is  no  more  than  about  25  per  cent  of  the  first-pass  absorption.  The 
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maximum  absorption  occurs  only  when  the  wall-to-wall  view  factor 
approaches  unity.  Subsequent  reflections  may  be  considered  but  they 
contribute  little;  the  second  reflection  amounts  only  to  about  6  per 
cent  of  the  absorption  by  direct  radiation.  Hence,  it  may  be  concluded 
that,  qualitatively,  a  maximum  error  due  to  imperfect  wall  emissivities 
can  be  estimated  if  the  first  few  reflection  processes  and  the  wall-to- 
wall  view  factors  are  considered. 

Imperfect  absorption  by  the  particles  is  more  complicated.  A 
first  approximation  for  absorptivities  less  than  unity  is  to  multiply 
the  ideal  result  by  the  absorptivity  of  the  particles.  The  results 
for  particle  absorptivities  less  than,  but  near,  unity  should  be  accur¬ 
ate  within  a  few  per  cent;  they  may  be  estimated  as  above.  Further  and 
more  precise  adjustments  are  complicated  by  scattering  patterns,  changes 
in  attenuation  coefficient,  and  view  factor  considerations,  to  mention 
some  of  the  more  important  variables.  Additional  study  is  recommended. 
The  problem  should  be  considered  along  with  back- radiation  to  which  it 
is  closely  related. 

Still  another  area  requiring  further  investigation  is  the  effect 
of  aerosol  temperature  variations  on  back- radiation  and  attenuation 
coefficients.  It  is  believed  that  this  work  can  readily  be  extended  to 
describe  these  more  involved  problems. 

Experimental  Results  and  Comparison  with  Theoretical  Development 

Previous  investigators  have  dealt  with  the  absorption  of  radiant 
energy  by  particles  clouds  and  have  presented  theoretical  formulations 
applying  to  their  problem.  Generally,  however,  the  calculated  results 

122 


served  as  corrections  to  other  problems  or  as  secondary  but  important 
factors  rather  than  as  a  study  of  primary  interest.  Reliance  was  made 
almost  completely  on  calculations  because  of  the  difficulties  in  ob¬ 
taining  experimental  data  for  radiation,  to  dispersed  systems.  This 
investigation  sought  to  measure  the  absorptances  of  selected,  well- 
characterized  aerosols  and  to  compare  the  results  with  theory  as  herein 
presented.  A  primary  effort  was  also  made  to  study  specifically  the 
effect  of  particle  size  on  the  absorptive  characteristics  of  the  aero¬ 
sols.  The  heat  transfer  data  obtained  have  already  been  presented 
graphically  in  Figures  26,  27,  and  28  and  tabulated  in  Appendix  B. 

The  results  as  shown  were  for  the  total  increase  in  heat  absorption 
due  to  the  addition  of  particles  to  the  aerosol  stream.  Further  con¬ 
siderations  were  required  before  the  heat  absorption  due  to  radiation 
could  be  evaluated. 

The  experimental  approach  was  predicated  on  the  enhancement  of 
radiative  heat  transfer  to  the  aerosol  system  and  the  reduction  to  a 
minimum  of  all  other  heat  transfer  processes.  Not  all  the  extraneous 
heat  transfer  effects  were  removed.  From  the  experimental  and  calcu¬ 
lated  data  of  Appendix  B  it  was  established  that  when  air  alone  was 
flowing  in  the  aerosol  system  an  average  of  7*8  Btu/min  was  absorbed. 

This  value  ranged  from  7.^  to  7.9  for  practically  all  tests.  This  heat 
transfer  was  sufficiently  small  for  radiation  to  contribute  significantly 
to  the  result.  Adjustments  in  the  final  result  were  required  to  account 
for  changes  in  convection  when  particles  were  added  to  the  system. 

Particles  in  a  flowing  gas  increase  the  convective  heat  transfer 
coefficient  above  that  for  the  gas  alone.  This  problem  has  been  studied 
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bably  reduced  slightly  by  the  increase  in  heat  transfer.  It  was 
assumed  than  that  the  mean  temperatures  were  very  nearly  equal  and  the 
heat  transfer  with  and  without  particles  was  essentially  determined  by 
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Figure  38.  Nusselt  Number  for  Particle-Gas  Suspensions  for  Experimental 
System. 
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(4A3) 


(Q  )  h 

=  _E 

(Q  )  h 

The  ratio  of  heat  transfer  coefficients  was  numerically  equal  to  the 

ratio  of  Nusselt  numbers  given  in  Figure  38.  Since  (q  )  was  found 

by  exqieriment  to  be  7.8  Btu/min,  (Q  )  could  be  estimated.  The  amount 

c  p 

of  energy  previously  calculated  by  equations  (3-6)  and  (3*15)  repre¬ 
sented  the  total  heat  absorption  due  to  the  addition  of  particles  to 
the  system  with  corrective  heat  transfer  included.  The  heat  absorbed 
by  radiation  alone  was  less  by  the  amount  (Q^)p  -  (\)a* 

The  data  for  total  heat  absorption  and  for  radiant  absorption 
alone  are  presented  in  Figures  39  and  4o.  The  results  were  replotted 
versus  k  rather  than  particle  cloud  concentration  because  the  latter 
presentation  was  more  general  and  was  useful  in  comparing  the  experi¬ 
mental  with  the  theoretical  results. 

Particle  Geometry,  Attenuation  Coefficient,  and  Associated  Quantities 
Of  basic  importance  were  particle  sizes  and  size  distributions. 
These  data  were  required  for  the  evaluation  of  such  quantities  as  the 
particle  surface  and  cross-sectional  areas,  S^,and  k,  respectively. 

An  attempt  was  made  to  derive  the  surface  areas  from  low  temperature 
gas  adsorption  isotherms  but  this  approach  was  found  to  require  much 
more  study  before  it  could  yield  reliable  results.  Surface  area 
determinations,  therefore,  were  calculated  from  direct  measurements  of 
statistical  particle  diameters.  These  size  distributions  are  given  in 
Appendix  A.  The  size  designations,  the  zinc  powder  excepted,  are  for 
the  screen  opening  size  employed  in  preparing  the  different  fractions; 
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Figiire  39*  Rate  of  Heat  Gain  by  Ferrous  Sulfide  and  Cupric  Oxide 
Aerosols  as  a  Function  of  the  Attenuation  Coefficient 
Before  and  After  Correction  for  Convective  Heat  Transfer 


Rate 


Figure  40.  Rate  of  Heat  Gain  by  Zinc  Aerosols  as  a  Function  of  the 
Attenuation  Coefficient  Before  and  After  Correction  for 
Convective  Heat  Transfer. 


the  quoted  diameters  are  not  the  particle  size  limits  for  each  powder. 
Imperfect  separation  and  particle  attrition  during  the  tests  produced 
broadened  distributions.  The  mean  particle  sizes  among  the  san^les 
were  significantly  different. 

The  geometric  quantity  most  significant  in  radiation  heat 
transfer  is  the  ratio  of  particle  cross  sectional  area  to  particle 
volume.  This  parameter  multiplied  by  the  particle  volume  concentra¬ 
tion  gave  the  attenuation  coefficient,  k.  The  calculated  results  for 
are  summarized  in  Table  2.  The  computational  procedure  is  described 
in  the  section  headed  Particle  Geometry  and  Measurement  of  Chapter  III. 

While  the  attenuation  coefficient  is  determined  solely  by  the  pro¬ 
duct  of  the  quantity  and  the  particle  volume  concentration,  two 
other  parameters  influenced  its  value.  One  was  the  decrease  in  concen¬ 
tration  with  increasing  aerosol  temperature;  the  other  was  the  increase 
in  hold-up  concentration  because  of  velocity  slip.  The  ten5)erature  rise 
of  the  aerosols  studied  decreased  the  calculated  flow  concentration  by 
about  10  per  cent.  The  flow  concentration  as  determined  experimentally, 
however,  was  not  a  measure  of  the  number  (or  volume)  of  particles 
actually  within  a  given  length  of  aerosol  conduit,  it  is  the  hold-up 
concentration  that  is  actually  needed.  Direct  experimental  data  were 
not  available  for  the  particle  sizes  employed  but  from  other  studies 
with  larger  particles^^  it  appeared  that  the  hold-up  concentration  was 
about  10  per  cent  larger  than  the  flow  concentration.  Since  precise 
data  were  not  available  and  the  two  effects  were  opposite  and  probably 
very  nearly  equal,  these  two  variables  were  neglected. 
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Table  2. 

Surface  Property,  S^,  Calculated 
Particle  Size  Distributions 

from 

Material 

Size  Designation 

S  X  10"^ 

V 

FeS 

53-88,  micron 

(i/ft) 

6.6io 

44-53 

10.55 

30-44 

13.75 

20-30 

15.08 

CuO 

53-88 

4.670 

44-53 

6.130 

30-44 

10.67 

Zn 

20-30 

18.03 

653.0 


1-10 


Radiation  Properties  of  Quartz  Enclosure  and  Particles 


The  quartz  aerosol  conduit  could  he  regarded  as  behaving  like  a 
black  enclosure,  although  the  radiation  passing  through  it  was  not 
black.  Radiation  within  the  furnace  but  outside  the  quartz  tubes 
approached  that  of  a  perfect  radiator  because  of  the  high  emissivity 
of  the  furnace  walls  and  because. the  furnace  was  nearly  a  complete 
enclosure.  This  energy,  however,  upon  passing  through  the  quartz  tubes 
was  screened,  and,  from  the  transmissive  properties  of  the  quartz  given 
in  Figure  20,  the  radiation  reaching  the  center  conduit  appeared  to 
have  been  principally  between  one  and  three  microns  in  wave  length  for 
a  black  body  furnace  temperature  of  2200^  R.  From  the  standpoint  of 
energy  distribution  with  respect  to  wave  length,  the  radiation,  then, 
was  far  from  being  black.  The  quartz,  however,  was  over  90  per  cent 
transparent  to  the  radiation  (one  to  three  microns)  reaching  the  inner 
concuit.  Radiation  crossing  the  tube  without  being  absorbed  continued 
through  the  wall  and  the  overall  effect  was  the  same  as  if  it  were 
absorbed  there  by  a  black  surface.  The  quartz  then,  because  of  its 
excellent  transmissive  properties,  behaved  as  though  it  were  black  in 
that  reflections  from  the  wall  back  to  the  aerosol  could  be  neglected. 

The  ferrous  sulfide  and  cupric  oxide  aerosol  particles  were  con¬ 
sidered  to  be  very  nearly  completely  absorbing.  They  were  black  in 
appearance  which  meant  that  radiation  in  the  visible  range  (0.3  to  0,8 
micron)  was  largely  absorbed.  While  this  degree  of  absorption  did  not 
necessarily  hold  for  longer  wave  length  it  is  probable  that  a  high 
degree  of  absorption  prevailed  up  to  three  microns  which  is  not  very 
far  in  the  infrared  region;  for  most  systems  the  infrared  region  must 
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be  considered  up  to  15  microns.  There  were  no  experimental  data  for 

ferrous  sulfide  to  indicate  the  absorptivity  of  this  mat^erial;  a  value 

51 

of  approximately  0.8  was  reported  for  cupric  oxide  for  very  low-tem¬ 
perature  radiation.  This  experimental  number,  however,  does  not 
correspond  very  well  with  the  short  wave  lengths  transversing  the 
aerosol  system. 

The  zinc  particles  had  an  absorptivity  of  approximately  0.045. 

Measured  values  of  0.04  to  0.05  tiave  been  reported  for  pure  zinc  sur- 

o  52 

faces  and  a  value  of  0.11  for  surfaces  oxidized  at  750  F.  The  zinc 
particles  used  in  the  experimental  tests  were  smoothed  and  polished  by 
extended  circulation  in  the  apparatus  and  should  have  had  an  emissivity 
very  near  the  minimum  for  the  pure  material. 

Comparison  of  Theoretical  and  Experimental  Findings 

The  ferrous  sulfide  and  cupric  oxide  aerosols  were  similar  in 
particle  size  ranges,  concentrations,  and  radiation  properties.  Their 
radiation  absorption  characteristics  were  likewise  similar.  A  joint 
comparison,  therefore,  was  made  of  the  experimental  results  and  the 
theoretical  calculations.  The  correlation  is  given  in  Figure  4l.  The 
experimental  data  are  corrected  for  convective  heat  transfer  as  pre¬ 
viously  described;  without  the  adjustment  the  experimental  results 
averaged  about  twenty  per  cent  higher  than  the  calculated  values  which 
were  derived  for  totally  absorbing  particles.  The  absorption  data  are 
given  as  the  fraction  of  energy  absorbed  over  energy  emitted  since 
this  ratio  (or  efficiency)  is  of  more  practical  importance  than  merely 
the  heat  absorbed  by  the  particular  aerosols.  The  abscissa,  k,  of  the 
plot  also  has  general  significance  in  that  its  value  is  jointly  deter- 
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Figiire  Ul.  Theoretical  Curve  and  Experimental  Data  for  Ferrous  Sulfide 
and  Cupric  Oxide  Aerosols  Relating  Absorption  Efficiency 
and  Attenuation  Coefficient. 


133 


mined  by  concentration  and  the  parameter,  S^,  which  in  turn  is  a  func¬ 
tion  of  particle  size  and  size  distribution. 


The  correlation  between  the  experimental  and  the  calculated  data 
appears  good.  The  data  points  for  the  two  materials,  however,  group 
above  the  theoretical  curve  for  the  ferrous  sulfide  and  below  it  for 
the  cupric  oxide.  This  grouping  may  be  due  to  different  convective 
heat  transfer  characteristics  for  the  two  aerosols  because  of  particle 
shape.  The  ferrous  sulfide  particles  were  irregular  and  the  cupric  oxide 
particles  were  more  nearly  spherical  or  similar  in  shape  to  smoothed 
pebbles  in  a  creek  bed.  It  is  believed  that  the  experimental  data 
establish  that  particle  cloud  radiation  properties  can  be  properly 
evaluated  provided  the  particle  emissivities  are  high,  the  particle 
size  distributions  are  accurately  determined,  and  scattering  effects 
are  insignificant.  These  conditions  have  been  assumed  in  previous 
studies  but  were  really  never  tested  by  direct  experimentation. 

The  results  should  be  applicable  to  processes  where  the  particle 
dimension  is  equal  to  or  greater  than  the  wave  length  of  the  radiation 
impinging  on  the  particle.  For  black  body  radiation  from  a  2000°  F 
source,  the  particles  shovild  be  around  15  microns  and  greater,  the 
particle  surfaces  should  be  highly  absorbent .  If  the  radiation  source 
has  a  black  body  temperature  of  10,000°  F,  or  more,  particle  sizes  down 
to  around  0.5  micron  may  be  reliably  treated. 

The  data  for  the  zinc  aerosols  were  compared  with  the  theoretical 
development  in  a  manner  similar  to  that  for  the  ferrous  sulfide  and 
cupric  oxide  aerosols.  There  was  the  difference,  however,  that  the 
theoretical  curve  had  to  be  corrected  for  the  zinc's  very  low  absorpti- 
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vity;  in  the  former  comparison  the  particles  were  apparently  very  nearly 
hlack.  While  there  are  many  factors  to  he  considered  for  reflective 
particle  surfaces ^  a  first  approximation  was  made  by  multiplying  the 
theoretical  result  by  the  apparent  emissivity  of  the  zinc  particles. 

The  resulting  curve  is  shown  in  Figure  42  with  the  experimental  data 
points.  The  approximation  seems  good  for  the  high  k  values.  The 
agreement  is  uncertain  for  more  dilute  suspensions  since  experimental 
points  were  not  obtained  below  k  =  59* 

General  use  of  the  theoretical  approximation  is  by  no  means 
recommended.  The  agreement  for  the  high  k  values  is  probably  due  to 
the  fact  that  absorption  took  place  principally  at  the  outer  region 
of  the  particle  cloud,  and  the  small  metal  spheres  produced  very 
pronounced  back-scattering.  Hence,  where  the  incoming  radiation  en¬ 
countered  the  particles  near  the  outer  limits  of  the  aerosol  it  must 
have  been  reflected  away  from  the  aerosol  without  significant  multiple 
particle  scattering  for  each  ray.  The  result  then  would  be  essentially 
that  of  a  smooth  zinc  metal  surface  absorbing  0.043  of  the  incoming 
radiation  and  reflecting  the  balance. 

While  the  estimated  theoretical  line  may  be  incorrect  by  a  large 
percentage,  its  use  for  a  high  concentration  of  metal  spheres  appears 
good.  At  lower  k  values,  or  more  properly  kR  values,  the  result  may 
be  used  for  estimation  purposes  where  no  other  more-refined  approach 
is  available.  This  should  be  the  subject  of  further  investigations  in 
the  future. 
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CHAPTER  V 


CONCLUSIOHS  AM)  RECOlVDyiENHATIONS 

The  following  statements  siimmarize  the  more  important  results 
of  this  investigation: 

1.  A  combined  analytical  and  ninnerical  procedure  was  developed 
to  solve  the  basic  radiation  equation  describing  the  radiant  heat  trans¬ 
fer  from  a  black  cylindrical  enclosure  radiating  uniformly  to  a  blacky 
evenly  dispersed  particle  cloud  within  from  which  back  thermal  radiation 
and  radiation  scattering  were  negligible. 

2.  A  generalized  graph  was  prepared  from  theoretical  concepts 

which  describes  the  rate  of  heat  absorption  by  the  particle  cloud  in 
terms  of  cylinder  wall  radiation  intensity^  cylinder  length  and  radius, 
and  particle  cloud  attenuation  constant.  The  results  were  expressed  as 
non-dimensional  quantities  with  ordinate,  L/R 

as  the  abscissa,  and  kR  as  a  parameter, 

3.  The  theoretical  results  for  a  given  kR  value  may  be  extra¬ 
polated  to  very  large  or  very  small  L/R  values  since  the  curves  in  these 
regions  become  very  nearly  linear, 

4.  Interpolation  between  kR  curves  may  be  accurately  accomplished 
by  replotting  the  data  on  a  log-log  grid  with  L/R  as  a  parameter.  The 
curves  so  obtained  are  smooth  and  for  the  most  part  are  linear,  except 
near  total  absorption. 

5.  For  L/R  ratios  greater  than  unity,  the  absorption  efficiency 


157 


of  a  particle  cloud  high  for  kR  values  of  unity  and  larger. 

6.  The  absorption  efficiency  decreases  rapidly  with  decreasing 
L/R  below  unity. 

7.  All  parametric  curves  approach  a  straight  line  with  a  slope 
of  45^  at  high  L/R  values  because^  for  these  conditions,  the  radiation 
system  is  behaving  as  though  it  were  an  infinitely  long  cylinder. 

8.  Equations  were  developed  expressing  the  rate  of  heat  absorp¬ 
tion  by  an  aerosol  while  in  the  entrance  and  exhaust  sections  adjoining 
a  heated  length. 

9.  An  experimental  method  was  developed  for  measuring  the  absorp 
tion  of  radiant  energy  by  particle  clouds  when  their  temperatures  are 
low  relative  to  that  of  the  source,  normally  convective  heat  transfer 
is  dominant  under  these  conditions  compared  with  radiative  heat  transfer 

10.  It  was  confirmed  experimentally  that  the  radiation  proper¬ 
ties  of  a  dispersion  may  be  characterized  by  its  particle  concentration, 
the  total  particle  cross  sectional  area,  and  the  total  particle  volume 
when  the  latter  are  determined  from  precise,  statistical  measurements 
from  particle  size  histograms.  Uncertainties  arise,  only  when  radiation 
scattering  is  significant  and  the  particles  are  not  black. 

11.  The  theoretical  development  should  be  applicable  in  deter¬ 
mining  the  dose  rate  for  an  aerosol  being  exposed  to  nuclear  radiation 
provided  the  appropriate  cross  sections  for  the  aerosol  particles  are 
known. 

While  the  results  of  this  investigation  contribute  to  an  under¬ 
standing  of  thermal  radiation  to  particle  clouds,  there  are  numerous 
other  problems  yet  to  be  solved.  Among  the  more  important  may  be 
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listed  the  following: 

1.  An  investigation  of  back  radiation  from  the  aerosol  to  the 
wall  should  be  made  using  the  approach  initiated  by  this  study.  The 
results  should  be  extended  to  include  radiation  scattering ^  and  the 
effect  of  non-black  particles. 

2.  The  effects  of  axial  temperatures  variations  should  be  in¬ 
vestigated  and  extended  to  include  radial  temperature  profiles, 

3.  A  calculation  of  the  effect  of  non-black  enclosures  on  the 
radiation  process  should  be  made. 

4.  The  experimental  technique,  perhaps  with  modifications,  should 
be  applied  to  systems  of  different  length-to-radius  ratios  to  test  the 
theoretical  work  more  thoroughly.  For  very  short  length-to-radius  ratios 
a  non-flow  system  with  a  radiation  section  constructed  similarly  to  the 
one  employed  in  this  study  would  probably  give  most  satisfactory  results. 

5.  Radiation  to  translucent  and  partially  transparent  particles 
has  been  studied  only  briefly  and  should  be  investigated  further. 
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APPENDIX-A 


PARTICLE  SIZE  DISTRIBUTIONS 


2  10  20  100 

Diameter,  microns 

Figure  4^+.  Particle  Size  Distributions  for  Ferrous  Sulfide 
Particles,  53-88p  Screen  Size  Fraction. 
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Per  Cent  by  Number  Less  Than  Indicated  Diameter 


Figure  45 .  Particle  Size  Distributions  for  Ferrous  Sulfide 
Particles,  44-53P  Screen  Size  Fraction. 
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Diameter,  microns 

Figure  k6.  Particle  Size  Distributions  for  Ferrous  Sulfide 
Particles,  SO-UUvi  Screen  Size  Fraction. 


Per  Cent  by  Number  Less  than  Indicated  Diameter 


100 


Diameter,  microns 


Figure  48.  Particle  Size  Distribution  for  Cupric  Oxide 
Particles,  53-88u  Screen  Size  Fraction. 
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Number  Less  than  Indicated  Diamete 


APPENDIX-B 
EXPERIMENTAL  AND 


CALCULATED  HEAT  TRANSFER  DATA 
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Experiment  conditions:  Heat  exchanger  coolant^  water^  0.095  CFM,  and  0.195  CEM  for  last  tw  tests; 
air  flow  rate  for  aerosol,  3*22  CFM  at  1  atm,  70  rj  Indicated  furnace  temperature  2280  F. 

Total  did  not  include  heat  calculated  for  coolant;  see  text. _ 
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Table  4.  Experimental  and  Calculated  Data  for  53-88p.  Ferrous  Sulfide’ 

(continued) 
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Experiment  Conditions:  Heat  exhchanger  coolant^  air  6.83  CEM  at  70  E_,  1  atmj  air  flow  rate 
aerosol  2.80  CEM  at  70 ^E^  1  atm;  furnace  temperature,  2280  E. 

Energy  balance  between  inlet  and  heat  exchanger  exit* _ _ 
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